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The EDSAC (Electronic Delay Storage 
Automatic Calculator) 


1. Introduction. The EDSAC is a serial electronic calculating machine 
now in operation in the University Mathematical Laboratory, Cambridge, 
England. It works in the scale of two and uses ultrasonic tanks for storage. 
The main store consists of 32 tanks each of which is about 5 ft. long and 
holds 32 numbers of 17 binary digits, one being a sign digit. This gives 
1,024 storage locations in all. It is possible to combine two adjacent storage 
locations so as to accommodate a number with 35 binary digits (including 
a sign digit); thus at any time the store may contain a mixture of long and 
short numbers. Short tanks which can hold one number only are used for the 
accumulator and multiplier registers in the arithmetical unit, and for con- 
trol purposes in various parts of the machine. 

A single address code is used in the EDSAC, orders being of the same 
length as short numbers. Most orders consist of a functional part which 
defines the operation and a numerical part which defines a storage location; 
some orders, however, consist of a functional part only. The complete order 
code is as follows. 


EDSAC Order Code 


Order Explanation 

An Add the number in storage location m into the accumulator. 

Sn Subtract the number in storage location m from the accumulator. 

Hn Transfer the number in storage location m into the multiplier register. 

Van Multiply the number in storage location m by the number in the multiplier register 
and add into the accumulator. 

Nn Multiply the number in storage location m by the number in the multiplier 
register and subtract from the contents of the accumulator. 

Tn Transfer the contents of the accumulator to storage location m, and clear the 
accumulator. 

Un Transfer the contents of the accumulator to storage location m, and do not clear 
the accumulator. 

Cn Collate the number in storage location m with the number in the multiplier 
register, i.e., add a 1 into the accumulator in digital positions where both 
numbers have a 1 and a 0 in other digital positions. 

R2"* Shift the number in the accumulator places to the right, i.e., multiply it by 2. 

LZ 2" Shift the number in the accumulator n places to the left, i.e., multiply it by 2*. 

En If the number in the accumulator is greater than or equal to zero, execute next 
the order which stands in storage location m; otherwise, proceed serially. 

Gn If the number in the accumulator is less than zero, execute next the order which 
stands in storage location »; otherwise, proceed serially. 

In Read the next row of holes on the tape, and place the resulting 5 digits in the 
least significant places of storage location n. 

On Print the character now set up on the teleprinter, and set up on the teleprinter 
the character represented by the five most significant digits in storage loca- 
tion m. 

Fn Place the five digits which represent the character next to be printed by the tele- 
printer in the five most significant places of storage location n. 

Y Round off the number in the accumulator to 34 binary digits. 

Zz Stop the machine, and ring the warning bell. 
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Ordinary 5-hole telegraphic punched tape is used for input. Each row of 
holes represents a 5-digit binary number and the basic input operation is 
to transfer this number to the store. Similarly, the output mechanism is a 
teleprinter, and the basic output operation is to transfer a 5-digit binary 
number to the printer, and to print the corresponding character. The tele- 
printer code is chosen so that binary numbers up to nine are printed as 
the corresponding figures and a similar code is used for input. This enables 
the operation of conversion to and from the decimal system to be pro- 
grammed as part of the calculation. 

The purpose of the F order is to enable the operation of the printer to 
be checked. Apart from this, no special checking facilities are provided in 
the EDSAC, and it is left to the programmer to incorporate in the program 
such checks as he considers necessary. 

2. The Control Sequence. When the machine is in operation, orders are 
executed automatically in the order in which they stand in the store. How- 
ever, when a conditional order (E or G) is encountered and the condition is 
satisfied, the next order to be executed is the one which stands in the storage 
location specified in the conditional order. 

Count is kept of the orders as they are executed by means of a short 
tank—known as the sequence control tank—which has associated with it 
an adding circuit through which unity is added to the number stored in the 
tank each time an order is executed. During a conditional order (when the 
condition is satisfied) unity is not added, but, instead, the number in the 
tank is replaced by the numerical part of the conditional order. 

The control sequence of the machine falls into two parts. In stage I, an 
order is transferred from the location in the store given by the number in 
the sequence control tank to a short tank known as the order tank. In 
stage II the order in the order tank is executed. 

The various constituent units of the machine—store, arithmetic unit, 
input unit, output unit, sequence control tank, and order tank—, are con- 
nected together through gates, so that the interconnections proper to each 
successive part of the control sequence can be made by opening the appro- 
priate gates. The gates are operated by waveforms supplied either by the 
main control unit or by a part of the machine called, for the purpose of 
this description, ‘‘order interpreter.” In the machine itself the “order 
interpreter” consists of a number of separate units. In the diagrams that 
follow, wires which carry pulses are shown as continuous lines, while those 
which carry control waveforms are shown as dotted lines. 

Figure 1 shows the state of the machine during stage I of the control 
sequence. Here it is assumed that the storage location specified by the num- 
ber in the sequence control tank is in tank 1 of the main store. The wire 
along which the order passes from this storage location to the order tank is 
shown by a heavy line in the diagram, and similarly the wires along which 
the control waveforms pass are shown as heavy dotted lines. The location 
section of the ‘‘order interpreter’’ receives the number in the sequence control 
tank and emits a waveform which opens the output gate of Tank 1. There are 
32 orders (or short numbers, or a mixture of orders, short numbers, and long 
numbers) circulating in the tank. The gate is opened as soon as the required 
number becomes available and closed again when it has passed through. In 
this respect it is treated differently from the other gates in the machine 
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which remain open or closed during the whole of stage I; The difference is 
indicated in the diagram by the different form of shading used. The 
switching and timing necessary for the operations detailed above are 
carried out within the “order interpreter.” 

Stage II is very similar to stage I except that the type of operation per- 
formed is not restricted to a simple transfer from the store to the order 
tank. During stage II, the order tank is connected to the “‘order interpreter” 
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in place of the sequence control tank, and the numerical part of the order is 
dealt with in the same way as the number in the sequence control tank during 
stage I. A new feature, however, is that the function section of the ‘order 
interpreter” comes into action and interprets the functional part of the 
order; it proceeds to emit waveforms which (a) set up connections between 
the units ready for the passage of numbers and (b) place the arithmetical 
unit (or input or output units, as the case may be) in a state of readiness 
for executing the order. Figure 2 shows the various gates set up for the 
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execution of an add order calling for the transfer of a number from tank 2 
of the main store to the arithmetical unit. One of the control wires passing 
from the ‘order interpreter” to the arithmetical unit is shown as a heavy 
(dotted) line to indicate that the arithmetical unit has been prepared for 
an addition. 

The above description has been written with the case of an arithmetical 
order involving the store in mind. In the case of other orders (e.g., left or 
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right shift, or a conditional order), stage II is simplified in that the location 
section of the ‘‘order interpreter” is not used. For example, during stage II 
for a conditional order, the order gate is connected to the sequence control 
tank through the gate shown in the diagram and the whole order transferred 
from one tank to the other. When an order is in the sequence control tank, 
only the numerical part is used, the function part being irrelevant. 

3. Initial Orders. From what has been said, and from an examination of 
the order code, it will be seen that the input mechanism is controlled by 
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program orders. Unless, therefore, there are some orders in the store at the 
beginning of the computation, nothing can be taken in through the input, 
and the machine cannot start. For this reason, there is a sequence of orders, 
known as initial orders, permanently wired onto a set of uniselectors (rotary 
telephone switches). These orders can be transferred to the store by pressing 
a button. 

There is considerable latitude in the choice of the initial orders, although 
once they have been wired onto the uniselectors, it is not easy to change 
them. The initial orders used in the EDSAC at present enable orders punched 
in the following form to be taken in from the tape. First a letter indicating 
the function is punched, then the numerical part of the order in decimal form, 
and finally the letter F or D indicating, respectively, that the order refers to 
a long or a short umber. If the order has no numerical part, it is punched 
simply as a letter followed by F. Under the control of the initial orders the 
machine converts the numerical part of the order to binary form and as- 
sembles the order with the function digits and the numerical digits in their 
correct relative positions. 


M. V. WILKES 


W. RENWICK 
The University Mathematical Laboratory 


Cambridge, England 


V. Wirkes, “The design of a practical niap-epeed computing machine. The 
EDSAC. ” R. Soc., London, Proc., v. 195, 1948, p. 274-279. 


M. V. WILEEs, “Programme design fora high-speed automatic calculating machine,” 
Jn. of Set Inst. and Phys. in Industry, v. 26, 1949, p. 217-220. 


Wirkes & W. RENWICK, “An ultrasonic memory unit for the EDSAC,” 
Electronic Engineering, v. 20, 1948, p. 208-213. 


Convergence Rates of Iterative Treatments 
of Partial Differential Equations 


1. Introduction. The development of high-speed digital computers' has 
made feasible the numerical solution by iterative methods of some partial 
differential equations. The convergence rates of several such iterative 
methods are estimated here. It is found that with the familiar elementary 
iterative methods some quite simple problems require prohibitive computa- 
tional labor. 

The iterative methods here considered are related to the various forms 
of the SouTHWELL “relaxation method’ ‘ in that they involve successively 
applied local corrections to improve an approximate solution. However, 
these iterative methods are routinized in conformity with the requirements 
of automatic computers while the relaxation method is flexible and depends 
in an essential way on the skill of its practitioners. 

2. Reduction to Finite Difference Form. The iterative methods of succes- 
sive approximation considered here are, like the relaxation method, not 
directly applicable to partial differential equations (and associated boundary 
conditions) but only to the finite difference approximations to them derived 
in the customary way.’ For example, the LAPLACE equations, A¢ = 0, 
applicable within a region ® in the x, y plane may be approximated by the 
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difference equation: 


(1) Lox, ¥) = o(x + h, y) + o(% — h, y) + ox, 9 + h) 

+ o(x,y — h) — 46(x, y) = 0 
applied at those points (x = jh, y = kh), of a rectangular lattice which lie 
within the region ®. We are here concerned only with the rate of convergence 
to the solution of this set of algebraic equations and not with closeness of 
that approximation to the solution of the differential equation. In the follow- 
ing we will use this difference equation, applied within a rectangular region, 
as an illustrative example for each of the iteration methods considered. We 
denote the value assigned to ¢(jh, kh) in the mth stage of iteration by ¢7.:. 
The limit approached with increasing ” we denote by ¢;,,. For definiteness 
we take & to be the rectangular region 


R = {j = 0,1,---,d, k =0,1,---,q}; 
on the boundary of this region we assign fixed values of ¢. Thus 
$3, = O50 = Die for j = O0or p, k = 0Oorg. 
(2) Loin = Oj-1,% + bitte + Oi0-1 + O5e41 — 4654 = 0 
for j=1,2,---,p—1 
k=1,2,-+-,qg—1. 
In all of the methods here considered ¢/, is a first (guessed) approximation 
to $;,x. Each of the succeeding approximations, ¢', ¢?, ---, is calculated on 
the basis of its predecessors (or immediate predecessor) by some process 
which guarantees the convergence of ¢" to ¢. 

3. The Richardson Method. Following SHoRTLEY, WELLER and FRIED‘ 
we term the most elementary of the familiar iterative procedures the 
RICHARDSON method.* In this method the correction process applied to each 
¢” consists of the addition of a positive multiple of L@* (for each interior 
point). Thus 

Ot = Oe + oldie = oe 
(3) + aL ¢j-ie + Of4i,2 + O5.2-1 + $3,241 — 463.2] (interior points) 
Oye = dhe = die (boundary points). 
The error at each stage we denote by ¢},z, 
jt = Oj.2 — 3% 
By substitution in (3) we obtain the error recurrence relation, 
jt = 2+ alex (interior points), 
since L¢;,, = 0; or more briefly 
(4) e"t! = (1+ aL)e" (interior points) 
= 0 (boundary points). 

The most familiar form of the Richardson method is that obtained by 

setting a = }. Then equation (3) reduces to 


o54e = Feo + Olas + OF2-1 + 97,241] (interior points) 
= bj, (boundary points). 


(5) 
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As shown below, this form is not only numerically more convenient, by 
reason of the disappearance of ¢3,, from the right member, but is also in one 
sense the most efficient form. 

To determine the convergence rate of this process we expand ¢° in the 
eigenfunctions of the operator L subject to the boundary condition of (4). 


These are evidently,® 
(7,2 . . . =1,2,---,p-—1 
(6) i = sin (xrj/p) sin (xsk/q), 


= my hk --5qg— 7. 

The corresponding eigenvalue of L we denote by Ly,« 

L§ = Sat Set Sete —4ee 
=[sin (xr(j—1)/p)+sin (xr(j+1)/p) —2 sin (xrj/p)] sin (xsk/g) 
+[sin (ws(k—1)/q)+sin (ws(k+1)/q)—2 sin (wsk/q)] sin (xr j/p) 


(7) =[2 cos (xr/p) —2] sin (#rj/p) sin (xsk/q) 
+[2 cos (xs/q)—2] sin (xrj/p) sin (xsk/q) 
. =([2 cos (xr/p)+2 cos (xs/q) —4 62" =Le, nf", 
ence 


Le,» = — 4[sin® (xr/2p) + sin? (xs/2q) ]. 
All of these eigenvalues are negative. The smallest and largest in magnitude 
belong tor = s = 1 and tor = p — 1, s = q — 1, respectively. We denote 
these by Lo and Ly. 
(8) Lo = — 4{sin? (x/2p) + sin® (x/2g)] = — w°(p* + ¢*) 
Ln = — 4[sin? (x(p — 1)/2p) 
+ sin? (x(q — 1)/2¢g) ] >= —84+ w(p?* + ¢”). 

The eigenfunctions of L are also eigenfunctions of the iteration operation, 
K = (1 + aL), corresponding to the eigenvalues 
(9) Ko, = 1+ aLe,». 


Each error-eigenfunction component is multiplied by its Ky. at each 
iteration. Thus if 
at = p as a,, Pe 


8 


then 
(10) ek = DL Gye) Ke, 2 


With increasing m those components, e“, for which 
(11) |Ko,»|< 1 


diminish in importance, the more rapidly the smaller is | Ky,,.)|. This itera- 
tion method is convergent (in general) only if (11) is satisfied for all r, s. 
The ultimate convergence rate is determined by the maximum magni- 
tude of K, 


(12) K* = | Ke.»|- 
Since by (9) Kg.) lies in the range 
(13) Kn =1+ aln < Koy £< 1+ alo = Kyo 
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the equality signs holding for (r,s) = (1,1) and (p — 1,q — 1), K is de- 
termined only by these extremes, 


K* = max{|Kn|, |Ko|}. 


As a@ increases from zero, Ko drops slowly from unity, K,, drops rapidly from 
unity. Thus K* = Ky > 0 so long as Ky > — K,,. For greater values of a, 
K* = — K,,, hence K* then rises with increasing a. The smallest K* (hence 
the most rapid convergence) occurs where 


Ko =1+ aly = — Kk, = — (1+ aL,), 
hence for 


(14) a = — 2/(Lo + Ln) = }. 


For this optimum a-value the error eigenfunctions of longest and shortest 
“wavelength,” «” and «-+*-), decay at the same rate 


(15) K* = 1 — sin? (w/2p) — sin® (w/2q) 
= $[cos (x/p) + cos (x/g)] = 1 — w°(p* + g*)/4 


(the short wavelength error alternating in sign) while other errors decay 
more rapidly. 

The numerical processes carried out in this iterative approach to the 
solution of a Laplace equation may also be regarded as the successive steps 
in the solution of the heat flow equation, 

te] , . 
. = Ad@ (interior) 
@=6 (boundary). 


The parameter, a, then plays the réle of 
At At At 


(Ax)? (Ay)? 

This method of solution of the heat flow equation is unstable’ for a appreci- 
ably greater than }, i.e., for a time interval appreciably greater than 
3(Ax)? = }(Ay)*. For smaller a-values the solution approaches asymptoti- 
cally a stationary form, hence one satisfying the Laplace equation. 

The optimum property of a = } in the Richardson treatment of the 
Laplace equation is not peculiar to the boundary conditions here considered. 
If (L + 4) has an eigenfunction, ¢jz, belonging to the eigenvalue (L’ + 4), 


(L’ + Adeje = €j-1,e Hejas.e + €j,2-1 + ©2421) 


then e”’ = (—1)*+*e’, if it is consistent with the boundary conditions, is an 


eigenfunction belonging to the eigenvalue (L” + 4) = — (L’ + 4), 
(L” + Adeje = — (L’ + Mee = efiie + eppse $+ eee + ejay 
L”" =8—-L’. 


If the boundary conditions permit this reversal, e’ — ¢’’, then for each eigen- 
value, L;, there also occurs the eigenvalue — 8 — L;. Then Z) + L, = — 8 
leading to optimum convergence for a = }. 

For other linear equations, 


Lo = 0, 
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where L is the linear difference operator corresponding to the linear differ- 
ential operator under consideration, the Richardson method, 


(16) ort = g" + alg", 


may be made convergent provided all of the eigenvalues of L are of the same 
sign (and not zero). 
The optimum convergence rate is given by 


a= — 2/(Lo + Ln) 


(17) K°=1~alyo 1 — 2%, 
where 7 denotes the ratio, 
(18) 7 Lo/(Lo + L,). 


Here again the Richardson method is formally equivalent to the solution 
of a partial differential equation in one more variable, 
dg 
ag ~ 1% 
the solution being carried to a sufficiently great ¢ to make the rate of change 
with ¢ negligible, hence also 


Lo = 0. 


4. The Liebmann Method. In the LieEBMANN method® a correction 
process like that of the Richardson method is applied to each of the lattice 
points in succession in a regular pattern. The ¢-value so corrected is used in 
all subsequent operations in that iteration step. It may thus be termed a 
“continuous substitution method.” In its simplest form the lattice is scanned 
in the same direction along successive rows. Thus, as applied to the Laplace 
equation and boundary conditions described above, the Liebmann iteration 
process may be written, 


One = bhe + oLeytie + iene + Opt + Peas — 402] 
(19) (interior points) 


= bie (boundary points). 
If a is again given the value } this expression becomes 
(20) oe = Beith e + hans + OF +ef.e4 J, 


thus L¢;,, is brought to zéro (momentarily) at each of the lattice points in 
succession. In this form the Liebmann procedure may be regarded as a very 
mechanical application of the relaxation method. 
The iteration equation for the error 
jt = je — dik 
is obtained by substitution in (19), 
oe = 2+ afetie + €j41,2 + et + €j.t+1 = 4¢5,x] 

(21) (interior points) 

=0 (boundary points). 
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This error interation process can be written briefly as 

(22) ent! = K(a)e", 

where K(a) is a linear operator depending on the parameter a, but now 
not simply related to the Laplace (difference) operator, L. We again examine 
the spectrum of eigenvalues of K(a) and regard the greatest magnitude of 
these eigenvalues as a measure of convergence rate. If K is an eigenvalue of 
(22) then its eigenfunction ¢;,, (we suppress indexing of the eigenfunctions 


and eigenvalues) must satisfy the following equation, obtained by substi- 
tuting (22) in (21) 


Kejn = 6,4 + al Kepe + j410 + Kejein + 6,041 — 462], 


Cok = Cnt = Cj,0 = Cj,e = 0. 


(23) 
We seek solutions of the form 


(24) 6. = A’sin (4rj/p)-B* sin (xsk/q) 
Substituting in (23) gives, 


(25) (K — 1+ 4a)A‘B* sin (xrj/p) sin (xsk/q) 
= B* sin (xsk/q)[KA* sin (xr(j — 1)/p) 
+ A sin (xr(j + 1)/p)] + aA? sin (4rj/p) 
< [KB sin (xs(k — 1)/q) + B* sin (xs(k + 1)/q)]. 


To prevent the appearance of terms in cos (xrj/p) and cos (wsk/q) in the 
right member of (25) we must require 


(26) A? = B= K. 


Since multiplication of A (or B) by —1 is equivalent to the replacement of r 
by p — r (or s by g — s) we may, without loss of generality, take A = B. 
Equation (25) then becomes 


(A? — 1+ 4a) = 2aA[Ccos (xr/p) + cos (xs/q) ] 
or 


(27) A? = 2atA + (4a —1) = 0; ¢t = [cos (xr/p) + cos (xs/q) ]. 


The appearance of two values (in general) of A for each (r, s) would seem 
to give more than the (p — 1)(q — 1) possible linearly independent error 
eigenfunctions; however, the replacement of (r,s) by (p — r, g — s) merely 
changes the sign of f, hence multiplies the two A-roots by —1. Thus (24), 
(26), and (27) define just the (p — 1)(q — 1) linearly independent eigen- 
functions required to form a complete set. 


If we again take a = } the roots of (27) are 
t 
A=0, 2 


The value A = 0 corresponds to the complete removal of an isolated 
error occurring at.j7 = 1 or k = 1. The other root gives 


(28) K = A? = }f = }[cos (ar/p) + cos (s/g) P. 
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The greatest eigenvalue occurs both for (r = s = 1) and for (r = p — 1, 
s = q — 1), and is, 


= 3[cos (x/p) + cos (x/q) P= [1 — w°(p* + g*)/4P. 


Thus in each iteration cycle the most resistant errors are reduced as much 
as in two cycles in the Richardson procedure. (This doubling of efficiency is 
qualitatively plausible since the numbers from which each ¢7%' is com- 
pounded are in the mean half as ‘‘old” in the Liebmann as in the Richardson 
procedure.) 

5. The Extrapolated Liebmann Method. Since the eigenvalues, K, are 
all positive it seems plausible that increasing a should improve the efficiency 
of this process. The procedure so obtained we may term the “extrapolated 
Liebmann method.”’ Examination of equation (26) shows that an improve- 
ment in convergence rate can be so obtained. For a > } the constant term 
is positive. We may therefore distinguish two ranges of a, hence of the 
middle term. Where 


(30) a <4a—1; f < (4a —1)/e? 

the two roots of (26) are complex conjugates having the magnitude 
(31) |A|= (4a—1)!; |K|=4e—1, 

which increases with increasing a. For 

(32) (4a — 1)/e®& <P <4; a<} 


the two roots for A are real, one of greater magnitude than (4a — 1)! 
(but <1), the other less. Differentiating (27) gives 


pronen. ° ° on im _ 4 
(33) d 4 2 ; At-—2<0; A-—ta= + [of — (4a — 1)}}, 





thus in the range (32) the greater of the magnitudes of the roots decreases 
with increasing a. The optimum value of a@ is therefore that for which the 
range (30) of complex roots just covers the spectrum of ¢-values. For the 
largest and smallest t-values 


(34) t = + [cos(x/p) + cos (x/g)] = + tmax; 
thus for a we use the smaller root of 

(35) Ofna — 4a +1 =0. 

Then for all of the error eigenfunctions 

(36) |K| = K* = 4a —1. 

For large p and g this may be approximated by 

(37) K* = 1 — V2e(p? + ¢)}. 


With this procedure and with the optimum values of a the number of 
iterations required to produce a substantial improvement in a trial solution 
increases about linearly with p and g rather than quadratically as it does for 
the Richardson method or the Liebmann method using a = }. 
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The usefulness of the extrapolated Liebmann method is limited by the 
difficulty of determining the optimum a-value for more complex problems 
than the example considered here. It seems likely, however, that in many 
similar problems a considerable improvement in convergence rate can be 
achieved by a suitable choice of a and that an approximate optimum a-value 
can be found empirically without great difficulty. 

The Liebmann and extrapolated Liebmann methods have an advantage 
—for some types of machine applications—over the Richardson and similar 
methods in that they require carrying as machine ‘‘memory”’ no more than 
one complete set of ¢;,.-values. However, for use with punched card machines 
or other calculating machines with severely limited internal (rapid-access) 
memory this advantage is offset by the difficulty of retaining the newly 
calculated values for use in the succeeding point and (more particularly) for 
the adjacent point in the succeeding row (or column). For such machines 
it is more convenient to use procedures in which a ¢-value when calculated 
may be stored in the “external memory” until needed in the next iteration 
cycle. 

6. The Second-Order Richardson Method. An improvement in the con- 
vergence rate of the Richardson method comparable to that achieved for 
the Liebmann method by extrapolation may be gained by retaining for use 
in the calculation of ¢77;' not only ¢},,x (j’, k’ running over points neighboring 
j, k) but also ¢33'. This modification is suggested by the formal equivalence 
between the Richardson method and the solution of the time-dependent 


<2 " . . . . 
equation, ed = Ld. The modified iteration method is equivalent to the 
solution of an equation of second order in ¢, 


Hp dg 

an + AG “ERAS 

and may therefore be termed the ‘“‘second-order Richardson method.” The 
iteration process may then be written as 


O37 = O3.2 + aLdj,e + B(G}.2 — $32) (interior points) 


(38) an Bes (boundary points). 


We again denoted by ¢;,, the solution of the equation L¢;,, = 0 satisfying 
the boundary condition incorporated in (38) and by é,x the difference, 
$3, — $;,x- The error then satisfies the induction equation (38) with a zero 
boundary condition. An error-eigenfunction, ¢},z, of L belonging to the eigen- 
value, Lw, 

(39) Lé,x = L,j,x (interior points) 


jr = 0 (boundary points), 
is then multiplied in each iteration cycle by a factor, K, such that, 


Khe = Ky, + aL,K,"™ + B(K," ie K,")e 


or (dropping the subscript, v) 
(40) | K?-(iteL+AK+ 6 =0., 
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If all of the eigenvalues, L,, are of the same sign (which for definiteness we 
take as negative) and are in the range 


(41) 0O>L2L2 Im, 
then @ and 6 can be chosen to minimize 
(42) K* = max|K,| < 1. 


(Since we are here using a second-order recurrence relation for ¢* the two 
roots of (40) correspond to two separate modes of decay for each eigen- 
function. Both decay rates must be considered in (42).) For positive 8 we 
again have a range of values of aL within which the roots of (40) are complex 
and of magnitude 6}. It may be shown, as before, that outside this range the 
greater |K| decreases with increasing 8. Thus a and # should be chosen to 
make the range of complex roots just fit the range (41). This requires 


1+ aL) + 6 = 26 
1+ aL, +6 = — 263, 


hence 
a(Lo + Lm) = — 2(1 + B) 
aLo sad REE (1 ey 6})?, 
or 
(43) n = Lo/(Lo + Ln) = 3(1 — B)?/(1 + 8B); 


6 — 2p(1 — 2n)* +1 =0. 

Choosing the smaller root gives for the decay factor, 

(44) K* = pt = 1 — 2(1 — 2m)(n — 9?) S 1 — 2h 
For the rectangular Laplace problem 


n= wr(p? + ¢*)/8, 
hence 


(45) K* =1 — r(p* + ¢)/V2. 


Thus a comparison of the convergence rate of the second-order Richardson 
method (45) with the extrapolated Liebmann method (37) again shows a 
factor of two advantage over the latter. 

It is not in general possible, to the writer’s present knowledge, to effect 
a comparable further improvement in the Richardson method by using a 
third or higher order induction process. 

7. Applications to Eigenvalue Problems. The methods here described 
may be used to determine the fundamental eigenfunction of an equation of 
the form 

(L+)o =0 


with homogeneous boundary conditions if the smallest value of \? is approxi- 
mately known. The operator, L, is replaced throughout by L + Ni, approx. and 
an initial trial solution chosen to approximate the fundamental eigenfunction. 
Then the higher eigenfunction components of the trial solution will decay 
during the iteration process while the fundamental remains approximately 
unchanged (depending on the accuracy Of Xo, spprox.)- Similarly a higher 
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eigenfunction can be calculated by replacing L by (L + 2}, approx.)”- It is to 
be expected that the iteration method so obtained would be very slowly 
convergent since 7 is the square of the ratio of the largest to the smallest 
value of |Z + Xi. approx. | (¢ ¥ 0). 

8. Time Estimates. In order to indicate the order of difficulty of typical 
problems to which these methods are applicable we consider two problems; 
the solution of the Laplace equation and of the biharmonic equation in a 
square region, 

p=qe=QN. 
The possibility, for some forms of the boundary conditions, of treating the 
biharmonic equation by factoring into a Laplace and a Poisson equation, 
will not be considered. 

We assume, for definiteness, that the iteration process is to be continued 
until all errors are reduced by a factor of 10-*. The number of cycles of 
iteration required will then be approximately, 


Nmax = 6/| log K*| 
while the number of ¢-values to be corrected in each iteration is (NV — 1)*. 


We denote by 7 the mean time required for each arithmetic operation. 
Then the time required for the evaluation of one ¢j,; is T where 


Method T 
Laplace-Richardson 4r 
Laplace-Liebmann (a = 4) 4r 
Laplace-Liebmann (optimum a) 7r 
Laplace-2nd order Richardson Or 
Biharmonic-Richardson 127 
Biharmonic-Liebmann (with extrapolation) 157 
Biharmonic-Liebmann (no extrapolation) 127 
Second-Order Richardson 177 


The total calculating time is given in the following table for a 10 X 10 and 
for a 20 X 20 lattice, and the asymptotic form for large N 


Approximate Total Calculating Time 


N = 10 N =20 N Large 
Laplace-Richardson 9-10'r — 1.6-10% 11 Nr 
Laplace-Liebmann (a = }) 5-104 8-105 5.6 Nr 
Laplace-Liebmann (optimum a) 1.2-104r 1.1-105r 15 N*r 
Laplace-Second-Order Richardson 3-104 3+ 105r 40 N*r 
Biharmonic-Richardson 1.1-10'r 8-108 14 N*r 
Biharmonic-Second-Order Richardson 4-105, 7-10%r 48 N‘r 
Biharmonic-Liebmann ? oe ? 


For the biharmonic equation the ratio, 7 = Lo/(Zo + Lm), has been approxi- 
mated by the square of the corresponding value for the Laplace equation. 

To convert these values to true time estimates we may, quite crudely, 
approximate + by 10-° days for electromechanical computers (e.g., punched 
card machines) and by 10-* days for entirely electronic computers. It is thus 
seen that with a fairly fine mesh the calculating time required with the 
slower machines is uncomfortably large for the Laplace equation and pro- 
hibitive for the biharmonic equation if the normal Richardson method is 
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used. Even with the faster machines the time required for the solution of a 
biharmonic equation by the methods considered here is uncomfortably large 
unless the second order Richardson (or probably also the extrapolated 
Liebmann) method is used. It is clear that for many problems of interest the 
simplest iterative procedures will prove impossibly tedious even with the 
fastest automatic computers. 

The apparent likelihood that the extrapolated Liebmann procedure 
would prove more rapidly convergent and more convenient for electronic 
computers than the second-order Richardson method would seem to justify 
an experimental study with such a computer. 

The writer is indebted to R. H. MacNgat and W. E. Mine for en- 
couragement and helpful discussion. 

STANLEY P. FRANKEL 
California Institute of Technology 
Pasadena, California 


1D. R. Hartree, “The ENIAC an electronic computing machine,” Nature, v. 158, 
1946, p. 500-506. [MTAC, v. 2, p. 222] 

IBM Automatic Contr olled Calculator, IBM gen a New York, 1945. 

Anon., “The UNIVAC,” Goctronce Industries, v. 2, 1948, p. 

H. H. Goxpstine & J. von NeuMANN, Planning and Codbee : of Problems for an Elec- 
tronic Instrument, Institute for Advanced Study, Princeton, 1947. [MTAC, v. a p. 54-56] 

2 R. V. SOUTHWELL, Relaxation Methods in En = ayy: Science, Oxford, 1 

- ‘ E. T. WaittaKer & G. Rosinson, The Calculus of Observations, no dy a ‘Glasgow, 
1932. 
4G. SHortLey, R. WELLER, & B. Friep, “Numerical solution of Laplace’s and Poisson's 
equations with applications to photoelasticity and torsion,” Ohio State University, Studies, 
Engineering Series, Bull. no. 107, 1942. 

5L. F. RicHarpson, “The approximate arithmetical solution by finite differences of 
physical problems involving differential equations, with an application to the stresses in a 
masonry dam,” R. Soc., London, Phil. Trans. s. A, v. 210, 1911, p. 307-357. ‘‘How to solve 
differential © en approximately by arithmetic, ” Math. Gazette, v. 12, p. 415-421, 1925. 

*R. Courant, “Uber partielle Differenzengleichungen, Congresso ‘Internazionale dei 
Matematici, Atti, Bologna, v. 3, 1930, p. 83-89. 

7 This can be shown by the’ method used by H. Lewy in “On the convergence of solu- 
tions of difference equations,” Studies and Essays Presented to R. Courant on his 60th Birth- 
day, New York, 1948, uP. 211-214. 

8H. LIEBMANN, ie angendherte Ermittelung harmonischer Funktionen und kon- 
former Abbildungen, ” Bayer. Akad. Wiss., math.-phys. Klasse, Sitz., 1918, p. 385-416 
Further references to the Liebmann method are cited in footnote 1 of Shortley, Weller & 
Fried. See also MTAC v. 3, p. 350, footnote 3. 


On a Definite Integral 


The function 
fle) = [> w+ 2)exp (— wu 


is tabulated by E. T. Goopwin and J. Staton! [MTAC, v. 3, p. 483] from 
series expansions and by numerical integration of the differential equation 
satisfied by the function. 

The integral can be evaluated explicitly in terms of two simple tabulated 
functions. The writer used the LAPLACE transform? to evaluate the integral. 
This method will now be given. 

We may write 


(1) I@® = in (u + #4) exp (— u*)du = im (1 + v)— exp (— #*)dv 
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and if we let 
Lt) = f” 10 exp (— pod 


denote the Laplace transform of I(#), we have 


2) LUO) = [+ #10 + vo = 4 + pL + np] 


= $(p) 
so that by the inverse transformation 
1 e+io 
(3) 10) =55 [~~ 9@) exp (pndp 
e* ct+io ia 1 ct+io — . 
(4) “El ze In (s — tds +35 z-te**(1 + 2)—"dz. 


The value of these integrals may be found by referring to tables of 
transforms.** The result is 


(5) I(t) = 4e* f- u—e"du + ate f : exp (u?)du 
so that* 
(6) f(x) = } exp (— x*)Ei(x*) + F(x), 
where 
Ei(y) = fi u—e“du and F(y) = exp (— y*) f . exp (u*)du 


are well-tabulated functions.' By expanding this result in an ascending series 
it can be identified term by term with equation (6) of Goodwin and Staton’s 
paper. 

I have communicated this result to Dr. Goodwin who informed me that 
although not aware of it he was within one step of deriving it from equation 
(4) of his paper. The following derivation is due to Dr. Goodwin: 

Equation (4) gives 

f' (x) + 2xf(x) = xt — x. 
Thus 


(7) d[exp (x*) f(x) ]/dx = mt exp (x*) — x exp (x’). 

If this could be integrated between suitable limits it would give f in terms 
of Ei(x*) and F(x) but, in fact, it is impossible to choose such real limits and 
it was at this point that Dr. Goodwin stopped. However, equation (7) can 
be written 


(8) ¢| exp (x?) f(x) — x f exp(u')au} = — x exp (x’). 
Taking the indefinite integral 


(9) exp (x*) f(x) — a f ; exp (u?)du = — f ; u- exp (u*)du + c 


1 2 
(10) m—> © f tetdt + c. 
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Comparison of the known behavior of f(x) as x — 0, namely 


f(x) ~ — Inx — fy 


with the fact that Ei(x*) ~ In x* + 7 shows that the constant in equation 
(10) can be taken as zero if the lower limit of the integral is taken as — . 
Equation (6) follows immediately. 

This corresponds to the use of a limit of © when (7) is integrated so that 
it is not surprising that Dr. Goodwin missed the result. 

I should like to express my thanks to Dr. Goodwin for permission to use 
his method in the latter part of this paper. 

R. H. RitcHie 

University of Kentucky 
Lexington, Kentucky 

1E. T. Goopwin & J. Staton, “Table of /5” e*du/(u + x),” Quart. Jn. Mech. Appl. 
Math. v. 1, 1948, p. 319-326. 
pe. 2H. T. Carstaw & & J. C. Jaeger, Operational Methods in Applied Mathematics, Oxford, 
y, k ‘ 046 Pires, Applied Mathematics for Engineers and Physicists, McGraw-Hill, New 

or 
> ee Modern Operational Mathematics in Engineering, McGraw-Hill, 

New York: 1944 

‘FMR Index, p. 189-191, 220, 231. 

*An ambiguity i in this result will be discussed in the next issue—Ep. 
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723(A-Z].—Harotp T. Davis & VERA FisHEerR, A Bibliography of Mathe- 
matical Tables. Copyright, 1949, by Harold T. Davis, xxii, 286, 27 leaves. 
For sale University Bookstore, Northwestern Univ., Evanston, IIl., bound 
$4.00. Mimeographed one side of each sheet. 20.9 X 27.5 cm. 
This copyright work has been issued in an edition of only 50 copies in 
order to test the possible demand for a revised and more complete volume 
of this kind. The title page states that it was “prepared under the direction 


‘of’ Professor Davis ‘“‘with the assistance of” Miss Fisher. We shall try, in 


a general way, to articulate clearly what is to be found in this tentative 
edition and to make some (among many possible) comments, which may 
serve as useful suggestions in any later revisions. 

The volume is divided into four parts: 


I. Introduction, leaves i—xxii. 
II. Bibliography of Mathematical Tables, leaves 1-196 + 27. 
III. Index of Tables by Classification of Functions, leaves 197-262. 
IV. Index of Tables by Names of Authors, leaves 263-286. 


We shall begin by first considering the 223 pages of Section II. There are 
here about 3,680 entries dated 1475-1948. Titles of almost all periodical 
articles are given, which is a highly praiseworthy feature, but there is a most 
undesirable total lack of uniformity of treatment of such titles, as well as of 
the abbreviated forms for titles of periodicals, and other details. In the case 
of pamphlets or books there is a similar lack; sometimes the total number of 
pages in the volume is given but more often not; names of places of publica- 
tions are omitted in a number of entries, and generally unintelligible Latin 
forms are also to be found; there is no uniformity in the forms of names of 
authors (at least six of which are incorrectly spelled), and titles such as 
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“Major General” and “Sir” are given in some cases but wholly lacking in 
dozens of others; dates of publication are absent in a number of cases; in 
at least one case the library and location of a publication is given, and in a 
few cases some inkling is noted as to the contents of entries, but the exact 
pages in books, where material of interest in the bibliography appears, are 
practically never given; the date of a third edition of an important work is 
given as if it were the first, and similarly for a third tirage; there are many 
titles with errors, some of the most fantastic character. There are at least 
three cases where two works of each of three authors are each ascribed to 
two different people. A number of dates of publications are incorrect. Wrong 
tabular references to certain authors are to be found in at least two cases. 

The arrangement of titles is alphabetical according to authors, except in 
such cases as reports of the Table Committee of the BAAS (where the au- 
thors of tables are almost completely known), Harvard University, Compu- 
tation Laboratory, NBSCL, etc. For each of 20 letters of the alphabet, one 
to three extra pages (27 in all), numbered xa, xb, xc, with additional titles, 
are added. 

The list gives evidence that the method of compilation included the 
copying (not always accurately) of practically every entry of published 
material in: 

(a) FMR, Index, 1946 (with the addition of titles in the case of periodical 

references) ; 

(b) D. H. LEnMER, Guide to Tables in the Theory of Numbers, 1941; 

(c) Scripta Mathematica, v. 2 (1934), p. 91-93, 297-299, 379-380; v. 3 
(1938), p. 97-98, 192-193, 282-283, 364-366; v. 4 (1936), p. 101-104, 
198-201, 294-295, 338-340. 

(d) J. HENDERSON, Bibliotheca Tabularum Mathematicarum, 1926; 

(e) J. W. L. GiatsHer, Report of the Committee on Mathematical Tables, 

1873; and some titles from: 

(f) R. C. ARcHIBALD, Mathematical Table Makers, 1948; 

(g) MTAC, nos. 1-26, 1943-1949; 

(h) A. DE Morgan, “Table,” The English Compaen, Arts and Science 
Sect., v. 7, 1861. 


These sources certainly account for more than 90% of the titles. The 
reviewer has verified that of the 3,680 entries there are 1,085 in the letters 
A-F and that of these 94% of the entries were in (a)—(h). In the 6% of new 
titles 36 referred to tables, 21 to theory, and 6 to calculating machines 
(including slide rules). In the whole book there are references to over 130 
authors of material on calculating machines, and 16 references to graphical 
aids. It will thus be seen that the title of Section II is a decided misnomer 
and that the title of the work under review is only partially descriptive. 

Because of the nature of the compilation there are scores of titles referring 
to material of no possible current use and almost wholly inaccessible to nearly 
every reader. There is no such exclusion as was practiced in the FMR, 
Index. In its present form Section II is wholly unreliable as a historical or 
bibliographical guide. An enormous amount of labor will still be necessary to 
bring this fundamental Section into a form having appeal to scholars. 

With Section II we next naturally associate the 24 pages of Section IV 
and the 12 pages of Subject Classification, p. xi-xxii, of I. In IV is an alpha- 
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betical list of authors followed by indications of the kind of tables, or other 
contents, in the entries of I1, according to the useful scheme of classification, 
which is a great elaboration of the classification given by Professor Davis in 
his Tables of the Higher Mathematical Functions, v. 1, 1933, p. 5-12. Ina 
general way the new classification was that adopted by the NRC Committee 
on Mathematical Tables and Other Aids to Computation. 

In Section III (64 p.), under about 300 subject headings are references by 
authors to the entries in IV. 

In practically no case is there any indication of the range of any table. 
No warning is anywhere given concerning wholly useless tables, because of 
their gross inaccuracies. 

If, in the whole work, 6% of the titles are new items for consideration 
by the authors of the possible forthcoming new edition of FMR, Index, a 
valuable contribution to scholarship will have been made. Furthermore, the 
ready subject index of Section III may, with the aid of other works, be 
frequently useful to specialists. Misspelled ‘“‘Napierian” on p. xi and 202 
should be wholly eliminated. 

In the Introduction is a table showing the “distribution by centuries of 
3,410 contributions to tables and table-making,”’ for the period 1500-1947. 
The difference between 3,410 and 3,680 appears to indicate the number of 
titles published in 1948. For the century 1500-1599 two titles are indicated, 
presumably for Rheticus and Otho (listed as a separate title and not as part 
of the work of Rheticus) ; Regiomontanus, 1475, is not counted. The numbers 
are indicated for the following quinquennial periods; for the century 1600- 
1700 the total number of titles is given as 70; 1701-1800, 84; 1801-1900, 985; 
1901-1947, 2,262. 

It has not seemed worth-while to give chapter and verse for every state- 
ment made in this review. In a work of this kind a “prepared-under-the- 
direction-of”’ prop is indeed a frail one whereon to lean. 

R. C. A. 


724[A, B, D, E, G, H, I).—Fritz Empe, Tables of Elementary Functions. 
Second Edition. With 83 Figures. Tafeln elementarer Funktionen., Leipzig, 
Teubner, 1948, xii, 181 p. 16.3 K 24 cm. Price 11.60 German marks, 
bound in boards and canvas. 

This work, first published in 1940, was a great elaboration of the 78 pages 
devoted to Elementary Functions in the third edition (1938) of JAHNKE & 
Empe’s Tables of Functions. Professor Emde has the following preface in the 
present edition: 

“The Second Edition of the Tables of Elementary Functions (an almost 
unchanged reprint of the First Edition of 1940) should have been issued 
in 1944. But after having been printed all copies were destroyed at the book 
binder’s by bombs and fire during the war. It is only now possible to reprint 
this edition from the same manuscript. Pretzfeld, January 1948” 

An offset reprint of the 1940 edition was made at Ann Arbor in 1945, and 
this was reviewed in MTAC, v. 1, p. 384-385. 

The only difference which the reviewer noticed, in comparing this edition 
with the work under review, was the substitution, on p. 10, of a fine relief 
(Fig. 3) of the function z*, z being complex (instead of z = x” on logarithmic 
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scale for —1 < y < + 1) supplementing the altitude chart, Fig. 79, p. 158 
[the function z* = (re‘”)*+iv], 
; R.A &i 


725(B, L].—H. Bremer, Terrestrial Radio Waves. Theory of Propagation. 
New York, Elsevier Publishing Co., 1949, x + 343 p. 17 X 24.5 cm. 
Price $5.50. 

Page 44 contains 3D values of $(3x)!, where x is either S + } or S + 3 
and S = 0(1)5. Page 45 contains 3D values of $(3x)!, where the x are the 
first six positive roots of J2;3(x) — J—2;3(x) = 0 or Ji/3(x) + J_1;3(x) = 0. 
Page 69 is a relief diagram of Hi,3(x + iy), for — 10 < x < 10 and — 5 
< » < 5. The branch cut is on the negative real axis. 

Chapter VI. ‘Numerical computations and results.’’ The results are 
mostly given in the shape of careful diagrams of which there is a large 
number. Other similar diagrams can be found throughout the book. 

A. E. 


726[(C, L].—HarvarD UNIVERSITY, COMPUTATION LABORATORY, Annals, v 
22. Tables of the Function = and of its First Eleven Derivatives. Cam- 
bridge, Mass., Harvard University Press, 1949, xviii, 241 p. 19.6 x 
26.6 cm. Price $8.00. 

This interesting volume was prepared as an aid to the application of 

FourRIER transforms. It gives to 9D the functions 


g(¢) = o sing 


and g™(¢), for 2 = 1(1)11 and ¢ = 0(1/360)20x — 2/360. 

The argument ¢ is expressed in degrees, the interval being 3°. The table 
is conveniently arranged so that at one opening there are 14 columns, the 
extreme ones giving the argument and the 12 others giving g™(¢) for 
nm = 0(1)11. The functions are very moderate in their behavior. In fact the 
bulk of the table gives only 8 significant figures. Beyond p. 38 all values are 
less than a tenth. However, the functions are still very much alive at 207. 

Since the tables give successive derivatives, interpolation is most easily 
performed by TayLor’s formula, the second-order formula being adequate 
for 8D work, linear interpolation giving accuracy to within 4 units of the 


’ 6-th decimal. 





The table was computed on Harvard’s Mark I from an 18D table of 
sin @, the successive g’s were computed from the recurrence formula 


g(¢) = [=e ein ¢ _ ne>(6)] / ¢. 


The whole computation required less than two weeks of machine time. 

An introduction (p. xiv-xviii) by R. C. SPENCER includes a few of the 
applications of the table to the theory of Fourier transforms, although no 
actual examples are worked out. Briefly, there are two main uses which 
facilitate the transformation of a function defined or given by a power 
series (really a polynomial of degree < 12) or a Fourier series (really a trigo- 
nometric polynomial of degree < 20). These methods apply to functions 
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F(x) of the real variable x, which vanish outside a finite interval, say 
— 1 <x <1. Then the Fourier transform G of F is given by 


G(u) = i F(x)e*™*dx. 


In case 
F(x) _ L A,x", 
n=0 
then 
G(u) = 2D Ani-*g(u), 
n=0 
and in case 
F(x) = >> A,e**, 
n=0 
then 


G(u) = 2 > Ang(u + nz). 


Thus, the table is used in the first case by rows and in the second case by its 
first column to obtain directly the values of the transformed functions. 
D. H. L. 


727[D, P].—JuLes Gaunin, L. HoupaiL_e, & A. BERNARD, Tables pour le 
Tracé des Courbes de Chemins de Fer, Routes & Canaux. Nouvelle édition 
revue et corrigée. Nouveau tirage. Premiere Partie: Tables Trigonométriques 
. . . Deuxitéme Partie: Recueil de Coordonnées . . . Paris, Dunod, 1948, 
xlvii, 181 p., xiv, 182 p. 13.7 K 21 cm. 

This is a new tirage of a very old book; it seems to be practically identical 
with the edition of 1922, when the chief author Gaunin was already dead; 
there was another tirage in 1925. According to the Catalogue of the Bibli- 
othéque Nationale there were a 1919 edition containing 426 p., a two-v. 
edition in 1911; and the second ed. in 1904, 2 parts in one v. The first edition 
of Tables Trigonométriques pour le tracé des Chemins de Fer . . . , was pub- 
lished in Paris, Dunod, 1862, xxxii, 181 p.; the second part, Recueil de 
Coordonnées, by Jules Gaunin, L. Houdaille, & A. Bernard, Paris, 1896, 
xxvi, 176 p. 

The first part gives a 6D table of the six trigonometric functions, versed 
sine, and versed cosine, for a = 0(30’")90°. There are also several other 
columns with values of 2a, 180° — 2a, 90° — a, 6D values of ra/180° and 
of (90° — a)/180°. The introductory pages deal with material of the tables, 
with trigonometry, and with some practical problems. 

The second part is devoted to tables for determining coordinates of points 
on circular arcs. The tables on p. 1-120 are for finding ordinates correspond- 
ing to abscissae measured on a tangent to the arc to be determined. The 
tables on p. 121-154 are for solving similar problems with reference to a 
chord of the arc; and the tables, p. 155-156, for points corresponding to a 
prolongation of the chord. 

Of the five remaining tables, A-E, in A, for R = 100(5)600(100)- 
3000(500)4000, 9D values are given for 27R; lengths of arcs corresponding to 
central angles, 1°, 1’, 1”, i.e., 2rR/360, 2xR/21600, 2xR/1296000; and angles 
in seconds corresponding to an arc of 1”; also the values of 1000/R. 
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In B, for R = 1", corresponding to 1(1)100, interpreted as central 
degrees, or minutes, or seconds, are given the corresponding lengths of arcs 
in meters. 

In C, for arc of .01(.01)1(1)10(10)100(100)1000 meters, and R = 100(50)- 
600(100)1000, 1.1(.1)3(.5)4(1)6, 10, are given the corresponding angles to 
the nearest thousandth of a second. 

Tables D and E are for conversion of degrees to grades, and for grades 
to degrees. 

We have elsewhere referred to such easement curves as the clothoid 
which are the ones of importance in modern times in laying out railways and 
other routes; see MTAC, v. 3, p. 146, 452. 

R. C. A. 


728[F].—H. CHaATLAND, “On the euclidean algorithm in quadratic number 
fields,” Amer. Math. Soc., Bull., v. 55, 1949, p. 548-553. 
For each prime p of the form 24 x + 1 less than 2% = 16384 with the 
exception of » = 73, 97, 193, 241, 313, 337, 457 and 601 there is given a 
representation 


(1) p = qum, + gems, 


where 41, g2, 1, m2 are quadratic nonresidues of », and where gq; and gq» are 
odd primes dividing gm, and qm: respectively to odd powers. No such 
representations are possible in the exceptional cases of p listed above. 

No column headings indicate which numbers in the table are the q’s. 
These seem to be the first and third numbers on the right of the equal sign, 
except in the case of = 409 and 577. 

Erpés & Ko! showed that in case the prime ~ has a representation (1) 
the field K(p*) has no euclidean algorithm. The reason for the upper limit 2“ 
is that according to a theorem of DAVENPoRT,? no field K(p!) is euclidean 
for p > 2%. 

oD. e. LL. 
seen P. ce & Cu. Ko, “Note on the euclidean algorithm,” London Math. Soc., Jn. v. 13, 
2H Davenront, “Indefinite binary quadratic forms and Euclid’s algorithm in real 
quadratic fields,’ to appear in London Math. Soc., Proc. 


729[F].—S. D. Cuowta & J. Topp, “The density of reducible integers,” 

Canadian Jn. Math., v. 1, 1949, 297-299. 

A definition of a reducible integer is given in RMT 734. An alternative 
definition is the following. An integer is reducible in case the prime factors 
of 1 + nm? are all less than 2m. The authors in collaboration with J. W. 
WRENCH have examined the first 5000 numbers and find that about 30 per 
cent of them are reducible. This conjecture is unproved. The number of 
reducible integers in each of the 50 centuries is tabulated. The table also 
gives the number B, of integers » whose prime factors are less than 2m}. 
These numbers are shown to have a density 


1 — In2 = .3068528. 


The table has been completely recalculated by one of the authors and was 
found to contain 25 errata. Since the table contains only 12 lines it is more 
economical to reproduce it below in corrected form than to single out its 
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errata. The number A, of reducible integers <m is given on the left, not on 
the right as stated in the paper. 


0 1000 2000 3000 4000 
n A, Bp A, By A, B, A, B, A, 8B, 
1-100 0 657 31 43 29 «43 35 41 28 42 
101-200 29 «50 27.43 30 «42 28 43 28 40 


201-300 27.47 3344 23° 42 24 43 28 «(41 
301-400 26 «645 27.4 32. 39 32 43 31040 
401-500 31045 3145 27 44 28 «(41 27.42 
501-600 29 «45 2344 32. (39 34 «4 37 DB 
601-700 29 «44 27 «40 26 «643 24 40 33 Al 
701-800 29 «44 35 43 3241 30 43 3s #9 
801-900 27 «44 27.= 4S 27.43 29 40 30 «43 
901-1000 2342 31-39 29 «42 20 «(41 38 41 
Totals 280 463 292 427 287 418 284 416 315 408 1458 2132 


D. H. L. 


730[F].—L. GamBELLI, “Sui caratteri di divisibilita con una tabella dei 

coefficienti di divisibilita di tutti i numeri da 2 a 101,” Period. Mat., s. 4, 

v. 27, 1949, p. 109-116. 

This note gives a table of the absolutely least remainder of 10* on division 
by m for all values of and each integer m < 101. In case this remainder is 
negative it is printed in boldface type. The purpose of the table is to give 
rules for finding the remainder on division of a given number N by m by 
replacing N by a linear combination of its decimal digits, the coefficients of 
this combination being tabulated. These coefficients are periodic, preceded 
by a preliminary nonperiodic part whenever m is not prime to 10. 

As an example of the use of the table, for m = 37 the table gives the 
period 

1,10, —11, 1,10, —11,---. 
The corresponding criterion for divisibility by 37 is as follows: Let the digits 
of N, beginning from the left, be do, d, dz, ---. Then N is divisible by 37 if 
and only if 
dy + 10d; — 11d, + ds + 10d, — 11d, + --- 
is divisible by 37. 
D. H. L. 


731(F].—H. Gupta, “On a conjecture of Miller,” Indian Math. Soc., Jn., 

v. 13, 1949, p. 85-90. 

The well-known function yu(m) of Méstus, which is +1 or —1 according 
as m is a product of an even or odd number of distinct primes and which is 
zero otherwise, plays an important role in the theory of distribution of 
primes. STIELTJES conjectured in 1885 that the function 

Mi(x) = & u(n) 
n<t 
satisfies 
| Mi(x)| < x! 


whose truth would imply the truth of the celebrated RIEMANN hypothesis. 
The conjecture mentioned in the title of the present paper involves the sum 
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function 
M,(x) = LD Mi(n) 
nsz 

and asserts that M(x) < 0 for x > 3. The paper verifies this conjecture for 
x < 20000. A table is given of — “M(x) for x = 25(25)20 000. There is also 
a table of — A(x) = — M(x)/x for x = [100(100)20 000; 5D]. 

The author conjectures that A(x)/logx is bounded. According to a 
heuristic argument of Brun, 


A(x) = — 24 12/x+--- 


The present table fails to support this result since A(x) < — 4 for x near 
18500. 


D. H. L. 


1Vicco Brun, “La somme des facteurs de Mdébius,”’ my 10. Skandinaviske Mate- 
matiker Kongres, Comptes Rendus, Copenhagen, 1947, p. 40-53 


732[F].—M. Petrovicu, ‘“Elementarna posmatrania o rasporedu omaniihk 
prostikh brojeva,’’ [Elementary observations on the distribution of small 
prime numbers], Srpske Akad. Nauka, Belgrade, Glas, v. 189, 1946, 
p. 5-45. 


Page 43 contains tables of the number of primes not exceeding x and 
the number of these which are of the forms 6m+ 1 and 6m —1 for 
x = 50(50)1000. These functions are compared with certain approximating 
functions. 


733(F].—H. Tietze, ‘Tafel der Primzahl-Zwillinge unter 300 000,” Akad. 
d. Wissen., Munich, math. nat.-Kl., Sitz., 1947, p. 57-72 (1949). 


If p and p+ 2 are both primes, then (p, » + 2) is called a prime pair or 
a pair of twin primes. The table in this paper lists the prime pairs <300 000 
by giving » + 2 in each case. This makes 2994 such pairs. 

The author is apparently unaware of three recent papers by SuTtTon,! 
CHERWELL? and SELMER & NesHEIM*® [MTAC, v. 2, p. 210, 342]. This 
latter table agrees as far as it goes (200 000) with the one under review. The 
discrepancy noted between Sutton and Selmer & Nesheim in MTAC, v. 2, 
p. 342, is thus to be blamed on Sutton, who omits two prime pairs between 
70 000 and 80 000, one between 90 000 and 100 000 and one between 120 000 
and 130 000. 

On p. 58 the author gives a rather incomplete list of errata in KRAITCHIK’sS 
list of primes to 300 000. Missing errata,’ however, do not happen to in- 
fluence the author’s table. 

D. H. L. 

1C. S. Sutton, “An investigation of the average distribution of twin prime numbers,” 
Jn. Math. Phys., v. 16, 1937, p. 1-42. 

2 Lorp CHERWELL, ‘ ‘Note on the distribution of the intervals between prime numbers,” 
Quart. Jn. Math., Oxford s., v. 17, 1946, 2. 

3E. S. SeLmer & G. NESHEIM, “Tafel der Zwillingsprimzahlen bis 200 000," K. Norske 
Videnskabers Selskab, Trondhjem, Forhandlinger, v. 15, 1942, p. 95-98. 

4M. KraltTcnrik, Recherches sur la Théorie des Nombres, v. 1, Paris, 1924. 

bd On page 156 of the reviewer’s Guide to Tables in the Theory of Numbers, ‘Washington, 


1941, is given what was hoped to be a complete list of errata in itchik's list of primes. 
Tietze, however, notes two additional errata, namely: 


for p=252141 read p = 252143 
for p=297671 read p=297971 [MTAC,v. 2, p. 313] 
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734[F].—J. Topp, ‘‘A problem on arctangent relations,” Amer. Math. 
Monthly, v. 56, 1949, p. 517-528. 
There are two tables giving the representation of arctangents of integers 
and rationals as linear combinations with integral coefficients of arctangents 
of integers. The numbers 


(S) 3, 7, 8, 13, 17, 18, 21, 30, --- 


have the property that, if m be one of them, arctan m can be expressed in 
terms of arctangents of integers <m and not in (S). Thus 


arctan 18 = 3 arctan 1 — 2 arctan 2 + arctan 5. 


The utility of such representations has been pointed out by J. C. P. MILLER 
[MTAC, v. 2, p. 62-63, 147-148] in connection with the preparation of a 
table of the Gamma function of a complex variable. The study of such 
relations goes back to Gauss. The numbers of the set (S) are called reducible. 

Table I (p. 525) lists all reducible integers not exceeding 342 and gives 
for each of these 100 numbers the expression of arctan m in terms of the 
arctangents of irreducible integers. There is also given an auxiliary integer 
Cn permitting an immediate passage from an arctangent to an arccotangent 
relation, where c, is the new coefficient of arccot 1. Thus for m = 18, 
C, = 1 and 

arccot 18 = arccot 1 — 2 arccot 2 + arccot 5. 


Two errata have been supplied by the author 
n=183 for c,=1 read c, =0 


for 4(1) + (3) read 7(1) — (2) 
307 for t, =2 read fC, = — 2. 


n 


Table II gives for each prime p < 409 of the form a? + 5? (i.e., for 2 and 
for all primes of the form 4% + 1) expansions of arccot (a/b), (a > 6) asa 
linear combination of arctangents of irreducible integers. In addition there 
are given the numbers c for converting to arccotangents. With each # is given 
also the least positive integer m, for which 1 + n?, = mp, the quotient m is 
given in terms of its prime factors. 

D. H. L. 


735[G].—ALBERT SapDE. Sur les Chevauchements des Permutations. Published 

by the author, Marseille, 1949. 8 p. 

The tables (for m < 12) are of various classifications of the permutations 
of distinct elements associated with the foldings of a linear strip of postage 
stamps which are characterized as being without ‘“chevauchements.” A 
chevauchement is an interlacing of number pairs a, a + 1, and , b + 1 fora, 
b both odd or both even. Also, for  < 8, the permutations of m distinct 
elements with one element fixed are tabled according to the number of such 
interlacings. 

J. RioRDAN 
Bell Telephone Laboratories 
New York, N. Y. 
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736[I].—_R. E. GrEENwoop. ‘‘Numerical integration for linear sums of 
exponential functions,” Annals Math. Stat., v. 20, 1949, p. 608-611. _ 
The author considers numerical integration over a finite range, using 

nm +1 evenly spaced abscissae. He introduces coefficients which would 

make the integration exact if the integrand were a linear combination of 
functions of the form exp (jx), where j runs from 0 to m (positive case), or 

from —m to m, 2m = n (symmetric case). For both cases, for m = 1(1)6, 

the coefficients are given to from 5D to 9D; about 45 values in all. He com- 

pares the results with the usual NeEwron-Cotes method, for m = 4, and 
finds that his coefficients ‘‘compare favorably’ with Newton-Cotes for the 
functions 1/(x + 3), exp (—x?), x exp (x), x*, and exp (2.2x). 

J. L. Hopcgs, Jr. 
Univ. of California 
Berkeley, Calif. 


737(I, L].—H. E. Satzer & Ruts ZucKER, “Table of the zeros and weight 
factors of the first fifteen Laguerre polynomials,’ Am. Math. Soc., Bull. 
v. 55, 1949, p. 1004-12. 


For the first fifteen LAGUERRE polynomials, the zeros x; are given to 
12 decimal places, the weight factors 
a” = [n!/L,! (x) 2 /x™ 


and values of a; exp {x;™} to 12 significant places. The chief use of the 
table is for numerical integration in the semi-infinite range by means of the 
mechanical quadrature formula 


JP pede = E as exp (21) fle.) 


which is exact when e*f(x) is a polynomial of degree not exceeding 2” — 1. 
A. E. 


738[K].—H. J. Gopwin, “On the estimation of dispersion by linear sys- 
tematic statistics,” Biometrika, v. 36, 1949, p. 92-100. 


Let x1 £ x2 < +++ < x, be the order statistics of a random sample of 2 
from a normal population with cumulative distribution function 


F(x) = (on) f exp (—42")dt, 
and let y; = %i41 — x; (¢ = 1,2, --+,m — 1). Table 2 gives E(y,) and E(y.y;) 
to 5 or more decimals for m = 2(1)10. The table of E(y.y;) supersedes an 


equivalent table to 2 decimals by Hastincs, MosTELLER, TUKEY, & WINSOR 
[RMT 740]. To calculate Table 2, ¥(i) and (i, 7) defined by 


vO = [" Pett — Fe)]ae 


Wid = [Pe ["O - Foy danas 
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were used; these are tabled to 10 decimals in Table 1, for ¢ = 1(1)5 and 
j = i(1)10 — 4. Suppose next that the differences y; are defined as above but 
with F(x) replaced by F,(x), where F,(é) = F(é/c). 

Table 3 gives to 5 or more decimals the coefficients a; which minimize 


n—1 
the variance of d = >) ax; subject to the condition that E(d) = ¢, for 
i=1 
nm = 2(1)10. Table 4 lists the efficiencies (calculated as inverse ratios of 
variances) of the unbiased estimates of ¢ formed from the following statistics: 
the “‘best’’ linear estimate d, the sample mean deviation from the mean, the 
sample mean deviation from the median, and x,44:1 — x: for k = 1(1)[4n], 
n = 2(1)10; the efficiencies are given to tenths or hundredths of a per cent. _ 
HENRY SCHEFFE 
Columbia University 
New York 27, N. Y. 


739[K].—E. J. GumBeEt, “Probability tables for the range,” Biometrika, 

v. 36, 1949, p. 142-148. 

The author is concerned with the asymptotic distribution of the range 
in a sample of m independent equally distributed random variables (that is, 
the difference between the largest and the smallest values in the sample). The 
desired distribution density is given by ¥(R) = 2e-*K,(2e-"2*), where 
K,(x) is the Bessel function. On p. 145 we find the values of ¥(R) and 


v(R) = i (x)dx. 


The range is —3.2(.1)10.6. The number of decimals drops from 8 at the tail 
ends (R < — 2.9and R > 9) to5 in the central part (about — 1.1 < R < 3). 
[See also MTAC, v. 4, p. 21.] . 


W. FELLER 
Cornell Univ. 


Ithaca, N. Y. 


740[K|.—Cecit HaAstInGs, Jr., FREDERICK MOsTELLER, J. W. Tukey & 

C. P. Winsor, “Low moments for small samples: a comparative study 

of order statistics,’’ Annals Math. Stat., v. 18, 1947, p. 413-426. 

This contains tables of the means, variances, and covariances of the 
order statistics in samples of size ¢ 10 from a normal universe, a rectangular 
universe, and a special universe with very long tails. For the normal universe 
and the special universe values computed from asymptotic formulas are also 
given. Means are given to 5D in all cases and are believed to be accurate to 
one unit in the fifth decimal. For the normal universe standard deviations 
are given to 5D with a maximum error of 2 or 3 units in the fifth decimal, 
while variances and covariances are given to 2D, with a possible error of one 
unit in the second decimal (except in one or two cases in which the error may 
be two units). Variances and covariances for the other two cases are given 
to 5D and are thought to be accurate to the places given. In addition the 
correlation coefficients in all cases are given to 2D. 

AE AE ot 
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741[K].—Pa.LmeEr O. JOHNSON, Statistical Methods in Research, xvi + 377 p. 

New York, Prentice Hall, 1949, 15 X 22.7 cm. Price $5.00. 

This book has five tables in an appendix. Tables I, II, III, and IV are 
standard tables of areas and percentiles for the normal, #, x?, and F dis- 
tributions. These tables occur in most modern applied statistics books 
or tables.* 

Table V contains first and fifth percentile values of the distribution of 
the NEYMAN-PEARSON L, statistic which is used to test the hypothesis that 
k normally distributed populations have equal variances. This test was de- 
signed for samples of sizes 11, m2, ---, a. These tables were constructed for 
the case m1 = mz = +--+ = m = n, and have as arguments k = 2(1)10(2)30, 
m = 2(1)10, 12, 15, 20, 30, 60, ©. They are reproduced from the Statistical 
Research Memoirs v. 1 edited by J. NEYMAN and E. S. PEARSON, and have 
not been widely reproduced elsewhere. These Memoirs are no longer easily 
available. 

Examples of the use of several other tables are included in the text. 
References are made to those tables, but they have not been printed in the 
text. Some would not be easily available to most readers of this book. For 
example page 102, reference 22, is to a table in Sankhya, v. 4, 1938. 

FRANK MASSEY 
Univ. of Oregon 
Eugene, Oregon 

1 For example, see: 

R. A. Fisher & F. YaTEs, Statistical Tables for Biological Agricultural and Medical Research, 
a 3rd ed., 1948 [MTAC, v. 3, p. 360-361] from which Tables II and III were 
produced. 


G. W. SNEDECOR, Statistical Methods A _ to Experiments in Agriculture and Biology, 
—_ Iowa, 1940, from which Table IV was reproduced. For errata see MTAC, v. 1, 
p. 85-86. 


742[K].—P. B. Patnaik, “‘The non-central x?- and F-distributions and their 
applications,” Biometrika, v. 36, 1949, p. 202-232. 
A noncentral chi-square variable, denoted by x”, with m degrees of free- 
dom (d.f.) and parameter x may be defined as follows: Let x; (i = 1, «++, ) 
be m independent standard normal deviates, let a; (¢ = 1, ---, ) be # con- 


stants, and define x” = } (x; + a,)*. The distribution of x” is known to 
i=1 


depend only on m and \ = 5 a,*. Noncentral F, denoted by F’, with »; and 
1 


v2 d.f. and parameter \ may be defined as a random variable distributed like 
(vex’")/(v1x?), where x” is a noncentral chi-square variable with », d.f. and 
parameter X, x? is a (central) chi-square variable with v2 d.f., and x” and x? 
are statistically independent. Various approximations to the x” and F’-distri- 
butions are considered. There are seven small tables, six of which compare 
these approximations with exact values, one of which (Table 6) gives the 
power of a chi-square test at the 5% significance level, calculated from one 
of the approximations to the x”-distribution. Tables similar to Table 6 but 
more extensive, for the 1% and 5% levels, and not based on approximations 
to the distributions, have been published recently by EVELYN Frx.' The 
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following approximations are considered for x” with n d.f. and parameter i: 
(i) x” is approximated by x?/K, where x? is a (central) chi-square variable 
with » d.f. (in general, fractional), and the constants » and K are determined 
by fitting the first two moments; (ii) normal approximation; (iii) and (iv) 
two different series for the cumulative distribution function of x” in each of 
which (i) contributes the leading term. F’ with »; and », d.f. and parameter 
\ is approximated as F/K, where F has the (central) F-distribution with » 
and y, d.f., K and » being fitted by the first two moments. The result is the 
same as though the above approximation (i) were used for the x” in the 
numerator of F’. 
HENRY SCHEFFE 


1 Evetyn Fix, “Tables of noncentral x*,” University of California, Publications in 
Statistics, v. 1, 1949, no. 2, p. 15-19. 


743(L].—W. R. Assort, “Evaluation of an integral of a Bessel function,” 
Jn. Math. Physics, v. 28,1949, p. 192-194. 


t 
The integral v(t, nm) =2n f u-'J2,(yu)du occurs in the theory of trans- 
0 


mission lines.! The expression of v as a finite combination of Bessel functions 
is known, but the two expansions obtained in this paper, 


2 n 
v=1——D (2k — 1)Ju«(y) 
vt emi 
and 
tt (—)9(2n — p - 1)! 
=i — 

e=1— 2 iin — p — DKw — py 
are better suited for numerical computation. The coefficients of the latter 
expansion are given numerically for » = 1(1)9. For m > 9 the first expansion 
is more useful. 





An-»(7) 


A. E. 


1 See for example, M. F. GARDNER & J. L. BARNES, Transients in Linear Systems, Studied 
by the Laplace Transformation, v. 1, New York, 1942, p. 310-317. 


744[L].—MiLton ABRAMOWITZ, “Asymptotic expansions of Coulomb wave 
functions,” Quart. Appl. Math., v. 7, 1949, p. 75-84. 
The author considers the differential equation 


d 
a + (1 — 2n/p)y = 0, 


and gives series-expansions of its solutions for different relative magnitudes 
of p and » (both real), which may be used for computational purposes. 
Starting with WHITTAKER’s integral representations of the confluent hyper- 
geometric function, he defines two solutions y; and y2 of the differential 
equation which are related to Whittaker’s confluent hypergeometric function 
by the relations: 
— 20(1 + inn 
2i(1 — in)y2 


Wg, -4(2tp) 
W_ ia, -4(—2%p). 
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He obtains asymptotic series for these solutions in the case that 7 is bounded 
and p— . These series are used to compute a table of the first three zeros 
of the function 


pPo(p, n) = e™%y1 + ¥2 


for » = 0(.5)3. Three other expansions for p®o(p, 7) are obtained, one for 
the case p < 2n, one for p > 2n, and one for 27 < 1 and p bounded. For 
bounded p and 7— ©, an approximation in terms of Bessel functions 
is given. 

Maria A. WEBER 
California Institute of Technology 
Pasadena, California 


745[L].—C. G. Darwin, “On Weber’s function,” Quart. Jn. Mech. Appl. 
Math., v. 2, 1949, p. 311-320. 


WEBER’S equation is the differential equation 


du 
qe + (+3 — i2*)u = 0. 


The case z = xVi,m = — 3 + ia (x, areal) arises in various wave problems, 
and is studied here. Instead of utilizing the theory of confluent hyper- 
geometric functions, of which WEBER’s function is a particular instance, the 
author develops his theory ab initio. The even and odd solutions of the equa- 
tion are unsuited because they are nearly proportional to each other for 
large x. The solutions selected for numerical tabulation now in progress at 
Scientific Computing Service, Ltd. (London) on behalf of the (British) 
National Physical Laboratory are 


ur = 2-*{(Gi/Gs) uo + (2G3/G1) *u1} 
ur = 2-3{ (G1/Gs3) 4 = (2G3/G:) *u1}, 


Gi =| + 4ia)|, Gs =|T(3 + 4ia)|, 


and w% and ; are the even and odd solutions normalized to values 1 and x 
for small x. For positive a, uy resembles an exponential function near the 
origin, and becomes oscillatory when |x| > 2a}; for negative a, the function 
is oscillatory throughout. Convergent power series for small x, and asymp- 
totic series for large x are given, both valid for fixed a. Other asymptotic 
expansions, valid for large a, are also developed. 

Weber’s equation, with real x and a, has also been investigated by 
Maaunus,! and by CHERRY.? 


where 


— = = 
1W. Macnus, Deutsche Math.-Ver., Jahrsb., v. 50, 1940, p. 140-161. 
?T. M. Cuerry, Edinburgh Math. Soc., Proc., s. 2, v. 8, 1948, p. 50-65. 


746[L]:—G. H. Goprrey, “Diffraction of light from sources of finite dimen- 
sions,” Australian Jn. Sci. Res., s. A, v. 1, 1948, p. 1-17. 


The paper contains tables concerned with the diffraction of light at 
rectangular and circular apertures. In particular Table 1 (p. 7) is a table of 


I(x) = [Si(2x) — x—' sin? x]/x 
for x = [0(.1)15(.5)34.5; 5D]. 
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There are also tables of the differences 


I(x + 2.49r) — I(x) for ax/e = — 1.2(.1)1.3 
I(x + 3.4%) — I(x) for x/e = — 1.7(.1)0 


to 5D. Table 4 (p. 13) gives 4D values of 


Tie) = f* (29/(xt*) ae 
for x = 0(.1)15, where 


HU) =< {2 + in sin (nein d — sao} 


is STRUVE’s function of order unity. 
Table 5 (p. 15) gives the difference 


Tix +84)—I(x) for «x =[—4.2(.1) — 2.8;4D]. 
D. H. L. 


747(L].—E. T. Goopwin & J. Staton, “‘Table of Jo(jonr),”” Quart. Jn. Mech. 

Appl. Math., v. 1, 1948, p. 220-224. — 

Five decimal values of Jo(joar) are tabulated for m = 1(1)10 and 
r = 0(.01)1; jo, being the mth positive zero of the BEssEL function of the first 
kind Jo. These values were calculated to seven decimal places by interpola- 
tion from the B. A. Math. Tables, v. vi, of Bessel functions; they were then 
differenced in the r direction on the National Accounting Machine as a 
check, and all doubtful roundings off were examined. The tables should be 
useful in the numerical discussion of potential problems with axial symmetry. 


A. E. 


748[L].—D. R. Hartree, ‘The tabulation of Bessel functions for large 

argument,”’ Cambridge Phil. Soc., Proc., v. 45, 1949, p. 554-557. 

The use of auxiliary functions in tabulating, to simplify interpolation, 
is well known. The author pleads for the use of auxiliary independent vari- 
ables for the same purpose. As an illustration, he shows that linear interpola- 
tion in a suitably constructed table of only 41 entries should be sufficient to 
give x4Jo(x) and x'Y,(x) from x = 5 to © with an uncertainty of one unit 
in the seventh decimal. The auxiliary variable is x~*, and the auxiliary func- 
tions are either P(x) and Q(x) defined by 


Jo + i¥o = (2/xx)"(P + iQ)et=, 
or R(x) and x(x) defined by 
Jo + Vo = (2/xx)tRe*@-*/4”), 
Tables are given for (2/x)§P and (2/z)#Q, and also for (2/x)*R and xp for 


ax? = 0(.01).05. 
A. E. 
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749[L].—HARVARD UNIVERSITY, COMPUTATION LABORATORY, Annals, v. 13: 
Tables of the Bessel Functions of the First Kind of Orders Sixty-four through 
Seventy-eight. Cambridge, Mass., Harvard University Press, 1949, x, 
566 p. 20 X 26.6 cm. Offset print. Price $8.00. 

A summary of reviews in MTAC of earlier published volumes of the 
Annals is given in RMT 711. The volume under review is the eleventh, in 
the monumental Harvard series of tables of BEssEL functions of the first 
kind, giving J,(x) for m = 64(1)78, x = [0(.01)99.99; 10D], but prior to 
x = 37.16 all values of J,(x) to 10D are zero. J7s(48.78) is the first significant 
value of this function. 

All values in this table have not been previously published. Two more 
volumes in preparation will deal with the next 22 orders, and give also the 
values of J,(100) for » = 0(1)100. 

R. C. A. 


750[L].—HaArvVARD UNIVERSITY, COMPUTATION LABORATORY, Annals, v. 21: 
Tables of the Generalized Exponential-Integral Functions. Cambridge, 
Mass., Harvard Univ. Press, 1949, xxv, 416 p., 20 X 26.6 cm., $8.00. 
The functions tabulated in this useful volume are 


E(a,x) = f (1 — e“)u-dt, Es(a, x) = r e~“u— sin u dt, 
0 0 
and ; 


Ec(a, x) = i (1 — e cos u)u—"dt, 


where u = (a? + £)}. The integrals 
E(a, x) = f e“u-'dt and f u- cos u dt 


can be evaluated in terms of elementary functions and the functions tabu- 
lated in this volume. All these functions are related to the generalized sine- 
and cosine-integral functions tables of which appeared in two earlier volumes 
of the same series! [MTAC, v. 3, p. 479-482]; they are also related (in 
special cases) to the exponential integral function by means of the relation 


E(0, x) + Ei(—x) = logx + 7, 


where log denotes the natural logarithm and y is EULER’s constant. These 
functions may also be regarded as incomplete (modified) BEssEL functions, 
for instance 


E(a, 2) -_ K,(a). 


Asymptotic representations (for large x) can be obtained by the usual 
method of integration by parts. 

Generalized exponential integral functions occur in the solution of the 
wave equation when line sources are immersed in a dissipative medium, and 
were encountered in particular in antenna theory.. Their tabulation was 
undertaken, on the Automatic Sequence Controlled Calculator, at the re- 
quest of Ronotp W. P. Kine. First the integrands were tabulated, and 
then the integrals computed by numerical integration, partly by Weddle’s 
rule, and partly by two rules based on fifth degree polynomials. An estimate 
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of the error in the integrands and of that due to quadrature, together with 
a study of the difference sheets, led to an estimate of 3 X 10-* for the 
maximum error in the integrals. The values, including differences; were 
rounded off to 6 decimals; an accuracy of 5.1 in the seventh decimal place 
is claimed for all numbers printed. 

The volume contains a Preface by Howarp AIKEN, an Introduction of 
four sections of which I “‘The generalized exponential-integral functions,” 
II “Computation of the tables,” III “Interpolation” were written by J. 
OrTEN Gapp and THEODORE SINGER who also collaborated in preparing the 
coding and control tapes and supervising the calculation, and IV ‘“‘Applica- 
tions” by R. W. P. King and C. T. Tat. There are 8 tables each for a different 
value of the interval 4. In each table both a and x run from 0 to 49h at steps 
of h. Each table consists of 50 pages of fifty lines each; each page contains, 
for a fixed value of a and for x running from 0 to 49h, 6 decimal values of 
E(a, x), Es(a, x), and Ec(a, x) together with first forward differences in both 
the a- and x-direction, except that differences are omitted when they can be 
found in the earlier parts of the volume. In the eight tables, h is respectively 
0.001, 0.002, 0.005, 0.01, 0.02, 0.05, 0.1, 0.2. 

All in all, this is a very useful and excellently arranged volume. 

A. E. 


1 HARVARD UNIVERSITY CoMPUTATION LABORATORY, Annals, v. 18, 19: Tables of the 
Generalized Sine- and Cosine- Integral Functions, Parts I and II, Cambridge, Mass., 1949. 


751[{L].—C. W. Jones, “On a solution of the laminar boundary-layer equa- 
tion near a position of separation,” Quart. Jn. Mech. Appl. Math., v. 1, 
1948, p. 385-407. 

This paper is mainly concerned with coordinating the investigations of 
the problem mentioned in the title by S. GOLDSTEIN on the one hand, and 
by D. R. HARTREE on the other hand. Since rigorous proofs seemed too 
difficult, the discussion is largely computational. 

Goldstein! expands the potential of the flow downstream in the form 


v = PLfoln) + Efiln) + Pfo(n) + --- J. 
The f, satisfy differential equations of the form 


fi" — 3b fe" + Gr + 2) f’ — (7 + 3)nf, = G,, 
where G, depends on fo, ---, f,-1, and suitable boundary conditions. 
fo = 0°/6, fi = ayn’, fe = aan? — a*n5/15. 


The differential equation for f,, for r 2 3, cannot be integrated explicitly, 
and so Jones uses for f; a step-by-step numerical integration, starting with 
the asymptotic solution for large 7. His Table 1 gives f,/a:’, and the deriva- 
tive of this function, for r = 0, 1, 2, 3 and m = 0(.1)4, to a varying number 
of decimal places. f, and fs contain an unknown constant, and in these cases 
only those parts of the functions, and of their derivatives, which are inde- 
pendent of that constant are tabulated, for the same 7, in Table 2. There is 
also a Table 3 giving > #f,’ as function of » for several values of a,t, but 
this table, as the author points out, is of less general application. 
A. E. 


1S. Goipstein, “On laminar boundary-value flow near a position of separation,” 
Quart. Jn. Mech. Appl. Math., v. 1, 1948, p. 43-69. 
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752[L].—GEORGETTE DE NOCKERE, Tables Numériques des Polynémes de 
Legendre P,,,o(cos 0) et des Fonctions Associées P,, (cos 0) ainsi que de leurs 
Intégrales P* jusqu’a n = 15 et 7] = 4, pour l' argument 6 (colatitude) variant 
de degré en degré. Tableaux des latitudes et longitudes divisionnaires et valeurs 
des multiplicateurs pour le calcul des coefficients du développement en série 
de polynémes de Laplace par la méthode des compartiments équivalents, 
d'une fonction de deux variables indépendantes. Acad. r. de Belgique, Ci. d. 
Sciences, Mémoires, v. 24, fasc. 4, publ. no. 1592, 1949, 166 p. 16.1 x 
25 cm. Price 150 Belgian francs. 
The first table, 5D, A (p. 12-43), isof P,, oand of Pio = S’ Pn, o(cos @)d cos @, 
for = 0(1)15, 6° = 0(1)90. Here Poo = 1; Pio = cos6@, Pio = }$ cos? 6; 
P20 = 3(3 cos 26 + 1), etc. 
The second main table, 5D, A, (p. 44-151), is of P,, ; and of 


Pi; = S P,; (cos 0)d cos @ 
for m = 1(1)15, 7 = 1(1)4, 6° = 0(1)90. Here 
Pi, = sin 6, Pi; = 3(sin 26 — 26) + 4/z; etc. 


The constants of integration are always chosen such that for @ = 90°, 
P33 = 0. Because of this condition the values for { P,, ;(cos @)d cos @ given 
by G. Prévost, Tables de Fonctions Sphériques et de leurs Intégrales. . . . 
Bordeaux and Paris, 1933, p. 153*-155*, have to be checked for appropriate 
constants of integration. In calculating the integrals, computations were 
made either directly, or with the aid of the following recurrence relations 
of LrénarpD:! 


(n+ 2)(n+1-— i f Passs(cos 6)d cos 6 
= (n — 1)(n + 4) f Prrs(c08 0d cos 0 
— (2n + 1) sin? 6P,, ;(cos 6). 
G- fr. 4+1(cos 6)d cos 0° 
= G+ Nt Dnt 1 — A f Pass(c0s 0d cos 0 
+ 27 sin 6P,, ;(cos @). 
The original contributions of Miss-Nockere consist in computing: 
(i) the values of P,, ;(cos 0) for m = 9(1)15, 7 = 1(1)4, by the formula, 
Pris, (cos 0) = (2n + 1)sin 6 Pa, j;-1(cos 0) + Pas, ;(cos 6) 
(ii) the values of P%,> for m = 0(1)15, by means of the classical formula, 


(2m + 1) f P(cos 0)d cos @ = P,,:(cos #) — Py_:(cos @) + const. 


(iii) the values of P,,;, m = 1(1)15, 7 = 1(1)4. 


The values of P,, ;(cos 0) for m = 0(1)8, 7 = 0(1)4, were adapted from 
the 7-10D tables of TaLLovist, (a) ‘“Tafeln der Kugelfunctionen P,(cos @),”’ 
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1905; (b) “‘Tafeln der abgeleiteten und zugeordneten Kugelfunctionen erster 
Art,” 1906, Finska Vetenskaps-Societeten, Acta, v. 33, nos. 4; 9. 

Since the Preface to the volume was written by Prévost, and dated 1938, 
it seems probable that these tables were computed some years ago. On the 
title page is the statement: “Impression décidée 4 mai 1948.” 

The final section (p. 153-166) is devoted to the subject matter referred 
to in the latter part of the title and is an elaboration of ideas developed in 
Prévost’s work mentioned above. 

az <. -»- 
1A. Lrénarp, “Formules de récurrence pour les intégrales des fonctions adjointes 


des polynémes de Legendre,”’ Acad. d. Sciences, Paris, Comptes Rendus, v. 196, 1933, 
p. 1773-1778 (there’s a slip in the page reference here). 


753[L].—J. PACHNER, “Pressure distribution in the acoustical field excited 
by a vibrating plate.’’ Acoust. Soc. Amer., Jn., v. 21, 1949, p. 617-625. 


Notations: 2, is the m-th root of the equation 
JI m(iz)Im' (2) — Jm(2)Jm'(iz) = 0, 
Ann = €mJm*(Znn)Jm?(Smnt)e2™, where é9 = 2 and ¢, = 1 for m = 1,2, ---, 
Bas - 28" Jn (8San) J mn—1(Smn) / San 
Cas = a Snn) | (Gund a~1(Sae))- 

Table I: zn, to 1 to 3 decimal places for m = 0 and m = 0(1)4,m = 1 
and » = 0(1)3, m = 2 and m = 0(1)3, m = 3 and m = 0(1)2, m = 4 and 
n = 0,1, m = 5(1)7 and = 0. 

Table II: Ao, to 4 significant figures for » = 0(1)3. 


Table III: Am, to 4 significant figures for m = 1 and m = 0(1)3, m = 2 
and m = 0,1, m = 3(1)5 and m = 0. 

Table IV: B,,, to 4 significant figures for m = 0, 1 and m = 0(1)3,m = 2 
and m = 0,1, m = 3(1)5 and 2 = 0. 

Table V: C,,. to 4 decimal places for m = 0, 1 and m = 0(1)3. 

Table VI: 2mC,,, to 4 decimal places for m = 0 and m = 0(1)3, m = 1 


and m = 0,1, m = 3(1)5 and m = 0. 
Table VII: Table of values where the partial acoustic pressure falls 
to zero. 
There are also various graphs representing numerical values. 
E. A. 


754(L].—Cuar.es H. Papas & Ronotp Kino, “Input impedance of wide- 

angle conical antennas fed by a coaxial line,’’ IRE, Proc., v. 37, no. 11, 

1949, p. 1269-1271. 

This paper introduces the auxiliary functions {,(x) = ga(x) + ib,(x), 
which are related to the spherical Hankel function of the second kind h,“ 
by the equation: 

bleh on hi” (x) 
" he? s(x) = nx hy (x) , 
and gives values for the real and imaginary parts of ¢,(x) in the form of 
curves covering the range 1 < x < 15 and for g; to gap and }; to bi7. 
- P. T. Nims 





Chrysler Corporation 
Detroit 31, Michigan 
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755[L].—M. Rotuman, “Table of fo? Io(x)dx for 0(0.1)20(1)25,” Quart. Jn. 

Mech. Applied Math., v. 2, 1949, p. 212-217. 

The function f(x) = Jo? Io(t)dt, where I(t) is the Bessel function of 
imaginary argument, is tabulated for x = [0(.1)20; 8S], and the function 
e~*f(x) is given to 9S for x = 15(1)25. In both tables modified second and 
fourth differences are given. 

These tables extend considerably those previously published [MTAC, 
v. 1, p. 250, v. 3, p. 308]. 

D. H. L. 


756(L].—V. V. SoLopovnikov, “‘O primenenii trapetsoidalnykh chastotnych 
kharakteristik k analizu kachestva sistem avtomaticheskogo reguliro- 
vaniia” [On the application of trapezoidal frequency characteristics to 
the analysis of the behavior of systems of automatic regulation ], Avto- 
matika i Telemekhantka, v. 10, 1949, p. 362-376. 


To aid in the treatment of the transform 


g(t) = (2/n) 4 x f(x) sin tx dx 
the author tabulates (p. 370-371) the function 
(2/m) {Si(kt) — (1 — k)“[Si(t) — Si(kt) + f° (cos t — cos kt) ]} 


for [k = 0(.05)1, ¢ = 0(.5)26; 3D]. This function, which is the transform of 
a standard “trapezoidal function,” is also graphed for k = .8, .9, 1. [cf. 
RMT 726]. 


D. H. L. 


757(L].—A. Uxricu, “Die ebene laminare Reibungschicht an einem Zylin- 
der.” Arch. d. Math., v. 2, 1949, p. 37-41. 


The functions appearing in the integration of the boundary layer equa- 
tions according to the method of BLasius' and Howartn’ satisfy certain 
nonlinear differential equations, of which the author gives explicitly those 
for the functions tabulated in his paper. 

Table 1. Tables of f; and fs with first and second derivatives to 4 decimal 
places for 7 = 0(.1)4. 

Table 2. Tables of gs and hs with first and second derivatives to 4 decimal 
places for » = 0(.2)4. 

Table 3. Tables of gz, hz, kz, Zo, hs, ko, jo, G9 with first and second deriva- 
tives to 3 decimal places for » = 0(.1)4. 

Tables of f; and f; have been first given by H1iEMENz.‘ Howarth® im- 
proved the values of f; and gave first approximations for gs and hs, which were 
later improved by FROssLING.? The functions with subscripts 7 and 9 are tabu- 
lated here for the first time. They were obtained by numerical integration, 
by Adams’ method, of the nonlinear ordinary differential equations which 
they satisfy. 

A. E. 
1H. Brastus, “Grenzschichten in Fliissigkeiten mit kleiner Reibung,” Zschr. f. Math 
u. a ay 56, 1908, p. 1. ia eae f 

2*N. FROssLING, ‘‘Verdunstung, Warmiibertragung und Geschwindigkeitsverteilung 
bei zweidimensionaler und rotationssymmetrischer laminaren Grenzstrémung,” Lund, 
Sweden, Univ., Acta, Afd. 2, Arsskrift, v. 36, 1940, p. 1-32. 
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3 L. Howarth, “On the calculation of steady flow in the boundary layer near the surface 
of a cylinder in a stream,”” ARC Report no. 1632, 1934. 

4K. Hremenz, “Grenzschichten an einem in einen gleichférmigen Fliissigkeitsstrom 
eingetauchten geraden Kreiszylinder,” Dinglers Polytechn. Jn., v. 326, 1911, p. 321. 


758[L, V].—G. N. Warp, ‘“‘The approximate external and internal flow past 
a quasi-cylindrical tube moving at supersonic speeds,”’ Quart. Jn. Mech. 
Appl. Math., v. 1, 1948, p. 225-245. 

In the course of the work indicated in the title, the following two 
functions are used:—W/(z) which is the inverse LAPLACE transform of 
1 — Ko(p)/Ki(p), and V(z) which is the inverse Laplace transform of 
— p-'Ki(p)/Ki'(p), where K,, is the modified BEssEL function of the third 
kind. Five decimal values of both functions, for z = 0(.2)10 have been com- 
puted by the Admiralty Computing Service of Great Britain and are repro- 
duced here together with values of W(z) for z = — 1.8(.2)—.2 computed 
by the author and Miss Roussak at Manchester, England. 

Cf. also British Admiralty Report SRE/ACS 89, 1945 [MTAC, v. 2, 
p. 294-295 ]. 

A. E. 


759[V].—ZpENEK Kopat, Tables of Supersonic Flow Around Cones of Large 
Yaw. Technical Report no. 5, Massachusetts Institute of Technology, 
1949. xviii + 125 p., 19.5 X 26.8 cm. 


The present volume is the third in the well-known series of computations 
of supersonic flow past cones.! The title of these tables is somewhat mis- 
leading. The title ‘Tables of Supersonic Flow Around Cones of Large Yaw”’ 
compared with the previous volume “On Supersonic Flow Past Slightly 
Yawing Cones” suggests two limiting cases, one for large and one for small 
angles of yaw. This, however, is not the case and the present tables deal 
with the second order approximation, the previous one with the first order 
approximation in a development in powers of the yaw angle e. Both tables 
are based on partly unpublished results of STONE. 

The radial and normal velocity components u, v and the pressure and 
density p and p are written in the general form 


f= fot eX fracosmd +e fan cos nd + say 


where ¢ denotes the angle of yaw, ¢ the circumferential angle. The axes are 
fixed with respect to the flow; the origin is at the apex of the cone. The 
circumferential velocity component w, i.e., the component in the direction 
of increasing ¢ is given in the form 


g=eD giasinnd+ & DY g2,sin mp + ---. 
n=0 n=0 


In these series fo; fi. and gin; fo, and ge, stand for the zero, first order and 
second order terms respectively. The cos and sin series are due respectively 
to the boundary conditions which also lead to the result that only f,,; and 
£:,1 differ from zero in the first order and feo, fe,2 and geo, g22 in the second 
order terms. 

Hence, the tables consist essentially of tabulated values for the feo, fe 2; 
Z2,9 and ge 2 corresponding to the five variables u, v, w, p, p. The velocities are 
referred to C, the maximum velocity obtainable by isentropic expansion. 
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Values are given for 7 cone opening angles in the range from 10° to 50° total 
included apex angle. An additional short survey of the main results is in- 
cluded. The representation, print, arrangement, etc., is similar to the previ- 
ous volumes and of very high quality. 

In general four significant figures are given in the numerical results. It 
is stated that this required an accuracy of six digits in the zero approxima- 
tion (TAYLOR-MACCOLL solution) and five digits in the first order (STONE) 
approximation. It is pointed out that this accuracy is greater than present- 
day experimental methods warrant but kept for future improvements in 
experimental technique. On this point the reviewer does not agree: The 
necessary accuracy of computations of perfect fluid solutions has a natural 
limit due to viscous effects. Viscosity is bound to have an effect for flow past 
yawing cones especially at larger angles of yaw. Hence, this point together 
with limitations due to possible convergence difficulties in the series in « seem 
to govern the desired accuracy rather than experimental errors. The reviewer 
feels that the accuracy in these tables is very probably higher than necessary. 
Still it may be argued that a standard solution like the cone case should be 
carried out with too great an accuracy just in order to enable an observer 
to evaluate viscous effects. The main justification of computing the cone 
flow in such great detail is, of course, the fact that here one standard solution 
is provided for comparison with more approximate theoretical methods on 
one hand and for experimental exploration of viscosity effects on the other. 

H. W. LIEPMANN 
California Institute of Technology 
Pasadena, California 


1See MTAC, v. 3, p. 37-40, 197-198. 


760[V].—_T. Y. Tuomas, “Calculation of the curvature of attached shock 
waves,” Jn. Math. Phys., v. 27, 1949, p. 279-297. 


In a previous paper! the author considered the well-known problem of 
supersonic flow past a pointed, curved, two-dimensional object to obtain 
an expression for the curvature of the attached shock-wave in terms of the 
curvature of the stream-lines immediately behind the shock. The present 
paper contains a numerical application of the previous result to obtain the 
ratio of curvature of shock-wave to curvature of stream-line at the vertex 
of the body. The expression for the curvature then involves only the stand- 
ard ratios characterizing an oblique shock—density on the two sides of the 
shock, etc. The set of tables and graphs expressing the results of the calcula- 
tion, therefore, are simply those for the much-discussed flow against a 
wedge with the one additional quantity, the curvature ratio, included. 
_ Consequently, it is surprising to be told in §3 of the paper that to the semi- 
vertex angle zero we must associate the shock inclined at 90° rather than 
at the well-known Mach angle arccsc M. One is presumably to infer therefore 
that the author considers either the larger shock-inclination branch, of the 
well-known double solution to the problem, to apply to physical reality—or 
that one has a discontinuous jump from one branch to the other in the 
neighborhood of vertex angle zero. (So far as the reviewer is aware, neither 
of these phenomena has been observed although attempts have been made 
at devising experimental conditions under which the former might occur.) 
The mathematical error in the paper arises from the author’s statement of 
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the compression condition as: 
(density behind shock)/(density before shock) > 1 rather than > 1. 


The worst physical error would seem to be that of publishing the discussion 
accompanying the tables without first looking at some of the experiments 
to which presumably they are meant to be applied. 

If one ignores the discussion of §3 of the paper, and substitutes for it the 
caution that for a given semivertex-angle one should choose the smaller of 
the two shock-angles, the tables and graphs of the paper should be useful 
for obtaining the initial curvature of the shock in the case of a curved, 
pointed, two-dimensional object such as-an air-foil. The tabulation is in the 
form of 13 tables, each for a single Mach number (M = 1.05, 1.08, 1.12, 
1.18, 1.25, 1.35, 1.47, 1.63, 1.83, 2.12, 2.56, 3.24, 4.45) with argument within 
the table being the shock-angle. There is included the additional useful table 
giving the relation between vertex-angle, shock-angle, and Mach number for 
the greatest vertex-angle with attached shock at given Mach number. The 
accompanying graphs exhibit curvature ratio and vertex-angle as functions 
of shock-angle. 


RicHArD N. THOMAS 
University of Utah 
Salt Lake City, Utah 


1T. Y. Tuomas, “On curved shock waves,”’ Jn. Math. Phys., v. 26, 1947, p. 62-68. 
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In this issue references have been made to Errata in RMT 729 (Chowla 
& Todd), 733 (Tietze), 734 (Todd), 741 (Johnson); UMT 93 (Fukamiya). 


168.—F. S. Carey, “Notes on the division of the circle,” Quart. Jn. Math., 
v. 26, 1893, p. 332-371. 
Table IV, giving the coefficient of the 6-nomial sextic, has the following 
errata. This list is the result of a recalculation of the table. 


p Coefficient of For Read 
61 x —27 27 
109 xt 39 135 
181 x 13565 1685 
193 x —5182 —5184 
x 1936 1744 
229 x° —2103 187 
241 x 594 580 
373 x 381 380 
397 = 4960 — 5040 
433 x — 130032 —1728 
457 x 3561 3461 
103 x° 1773 1373 
127 x —977 —972 
151 x® 6547 6543 
163 x 21323 5023 
223 x —3276 5644 
x® —71228 4592 
331 x® 84429 84427 


EmMMA LEHMER 
942 Hilldale Ave. 


Berkeley 8, Calif. 
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169.—NBSMTP, Tables of Spherical Bessel Functions, v. 1, New York, 
Columbia University Press, 1947 [MTAC, v. 2, p. 308-309]. 


In the table of (x/2x)*J;(x) at x = 7.45 
for .12340 32451 read .12342 32451. 


GERTRUDE BLANCH 
NBS Institute for Numerical Analysis 
Univ. of California, Los Angeles 


170.—NBSMTP, Tables of the Exponential Function e*, New York, 1st ed. 
1939, 2nd ed. 1947 [MTAC, v. 1, p. 438, v. 2, p. 314, v. 3, p. 173]. 
P. 188, at x = 1.8784, for x = 1.9884 read x = 1.8784; and for 
e* = 6.54302 76384 25706 read 6.54302 76384 25796; for x = 1.9885, 1.9886 
read x = 1.8785, 1.8786. 


PauL ARMER 
Rand Corporation 


Santa Monica, Calif. 


171.—L. Scuwarz, “Untersuchung einiger mit den Zylinderfunktionen 
nullter Ordnung verwandter Functionen’’, Luftfahrtforschung, v. 20, 
1943, p. 341-372. [Translated by J. Lotsor, Cornell Aeronautical 
Laboratory, 1946. ] 

The following errors were found by differencing the tables and conse- 
quently the corrected values can be in error by at most two units in the last 
decimal place. 


Function r x for read 
SJA(r,x) 0.2 1.58 0.9_5297 0.995297 
JTe(d,X) 1.0 1.74 0.884505. 0.884547 
N.(A,x) 0.1 0.88 1.784990 1.784986 
N.(A,x) 0.1 0.96 1.861997 1.861965 
N.(A,x) 0.1 1.38 1.357677 1.357670 
N.(A,x) 0.2 0.44 0.967949 0.966929 
N.(A,x) 0.2 0.46 0.995782 0.995784 
N.(A,x) 0.2 0.58 1.155928 1.153056 
N.(A,x) 0.2 0.30 0.097532 0.097509 
N.(A,x) 0.2 0.44 _ 0.184975 0.184746 
N,(A,x) 0.2 0.58 0.289383 0.288526 
N,(A,x) 0.2 0.60 0.304306 0.304406 
N-(A,x) 0.7 1.48 0.831510 0.831310 
N,(A,x) 0.7 0.06 0.003768 0.004336 
N.(A,x) 1.0 1.54 0.005461 0.005036 
C.(A,x) 0.0 0.10 0.069995 0.070995 
C.(A,x) 0.0 0.58 0.064846 0.064865 
J.(d,x) 0.2 41 1.27530258 1.57530258 
N.(A,x) 0.9 4.9 1.388321 1.388821 


A. H. ROSENTHAL 


NBS Institute for Numerical Analysis 
Univ. of California, Los Angeles 
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172.—G. N. Watson, “A Table of Ramanujan’s function 7(m),’’ London 
Math. Soc., Proc., s. 2, v. 51, 1949, p. 1-13 [MTAC, v. 3, p. 468]. 
P.12, m = 847 for 38152 read 58152. 

This typographical error was noted when Watson’s table was put on 
punch cards and submitted to a series of checks. These cards are available 
in the UMT Fie. 

D. H. L. 


UNPUBLISHED MATHEMATICAL TABLES 


EDITORIAL NOTE: Beginning with this volume we are starting a collection 
of unpublished mathematical tables to be known as the UMT Fine. Authors 
of tables which have no immediate prospect of publication are invited to 
submit copies for deposit in UMT Fixe. Description of such tables will 
appear in UMT and photostat or microfilm copies will be supplied at cost 
to any reader of MTAC. Address tables or correspondence to D. H. LEHMER, 
942 HILLDALE AVE., BERKELEY 8, CALIFORNIA. 


9O[F].—R. A. LIENARD, List of primes of the form k-10° + 1 and k-10’ + 1 
for k < 1000. Manuscript in the possession of the author and deposited 
in UMT FILe. 

There are listed 117 primes of the form k- 10° + 1 and 109 primes of the 

form k- 10? + 1. 


R. A. LIENARD 
95 Rue Béchevelin 


Lyon, France 


91[K].—BALListic RESEARCH LABORATORIES. Aberdeen Proving Ground, 
Md., Probability Integral of Extreme Deviation From Sample Mean. 


H,(x) = (n/2x(n — 1)) f exp (— #/2n(m — 1))H,-:(t)dt, (Hi(x) = 1.) 


For normally distributed and ranked variables u; < uz < +++ < Usp, 
H,(nk) gives the probability that the extreme deviation from the sample 
mean, i.e., %, — &@ or %@ — x will not exceed & times the population standard 
deviation for random samples of size n, for m = 1(1)25. Each function is 
tabulated until it becomes sensibly unity. The interval of tabulation for x 
is as follows: 

1 for m = 2(1)7, A for n = 12, 13, 
.2 for m = 8(1)11, .8 for m = 14(1)25. 


The accuracy is 7D. 
The work was done on the ENIAC under the direction of J. V. Hot- 
BERTON. 


92[K].—BaListic RESEARCH LABORATORIES, Aberdeen Proving Ground, 
Md., Binomial Probabilities. 





I,(c,n —c+ 1) = Bn = c+ 1) [a — x)"*~*dx. 
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The probability of c or more successes in m trials is given by 
I,(c, m — ¢ + 1), where p is the chance of a success in one trial. 


nm = 1(1)200, c = O(1)n, p = 0(.01)1; 7D 


The work was done on the IBM Relay Calculator under the direction of 
M. LorKn. 


93[L].—FAsANoRI FuKaAmiyaA, Mathematical Studies of the High Frequency, 

Photostat filed in the Japanese Section of the Board of Trade, Technical 

Information and Documents Unit. Reference BIOS/JAP/DOC/1550. 

This paper contains a certain amount of descriptive text followed, it is 
stated, by 5 appendices. 

Appendices I and II consist of tables and graphs of certain functions 
(which are not of any interest being purely “ad hoc’’) connected with radia- 
tion from a sectoral horn. 

Appendix IV is stated to be a table of ‘zero points of the derivative of 
the Associated Legendre function.” In actual fact Appendix IV consists of a 
series of graphs of some unknown function. 

Appendix V is the most interesting part of the paper. This consists of 
tables of J,,(z) as follows: 


1 z = [0(.1)3;7D], n 


n ’ 
11(1)20, s=[0(.1)6;7D], n 


n 


2(1)10, 2 = [0(.1)10.7; 7D], 
11(1)20,. z = [6.4(.1)10; 5D]. 


(no explanation is given for the absence of J,(z) for m = 11(1)20, 
2 = 6.1(.1)6.3.) 

The tables were compared as far as possible with WatTson’s Bessel 
Functions! and, apart from differences of 2 or less in the last place, the 
following discrepancies occur. 


Watson Fukamiya 
J;(2.3) .1799789 .1899789 
J3(4.3) 4333147 4333470 
J5(1.0) .0002498 .0002475 


It is further stated that tables of Fresnel Integrals for argument s 0-50 
have been calculated, but were “destroyed in the event.” 

Unfortunately the reproduction of the copy seen was poor, only the 
three pages of text were typewritten, the tables being written. In addition 
appendix III was missing and no date was given for the work, nor was there 
any indication of where Fukamiya worked. 

Li. G. CHAMBERS 


Royal Navai Scientific Service 
London, England 


A comparison with the forthcoming B. A., Mathematical Tables, v. X 
revealed many errors; these are often only a unit or two in the last figure, but 
many are gross errors. Some effect several values of m for the same x and 
are obviously computational. It is noteworthy that errors are fewer in the 
range covered by Watson’s table! and are almost all in the end figure, yet 
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agreement with Watson is not complete, and a value known to be in error 
in Watson is correct in the table under review. 
J. C. P. MILter 
23 Bedford Square 
London W.C. 1 


1G. N. WatSsON, Treatise on the Theory of Bessel Functions, Cambridge, 1922, second 
ed. 1944 [MTAC, v. 2, p. 49-5i.] 


94[V].—BALLIsTICc RESEARCH LABORATORIES, Aberdeen Proving Ground, 
Md. Supersonic Flow past Cone Cylinders. [See Note 113.] 


AUTOMATIC COMPUTING MACHINERY 


Edited by the Staff of the Machine Development Laboratory of the National Bureau 
of Standards. Correspondence regarding the Section should be directed to Dr. E. W. 
Cannon, 225 Far West Building, National Bureau of Standards, Washington 25, D. C. 


TECHNICAL DEVELOPMENTS 


Characteristics of the Institute for Numerical 
Analysis Computer 


In January 1949, members of the staff of the Institute for Numerical 
Analysis' began the development and construction of a high-speed electronic 
digital computer. As of December 1, 1949, the central computer was approxi- 
mately eighty per cent completed. The group responsible for the building of 
this machine is composed of, besides the author, three engineers, three junior 
engineers, and four technicians. In addition, one mathematician is assigned 
to the coding and programming of problems to be run on the machine. 

Information is stored and processed in the computer in units called 
words, a word consisting of 41 binary digits. This word length is determined 
by the number of words which can be stored in the computer’s high-speed 
memory. 

Words in the machine sense may represent (1) numerical information, 
(2) instructions to the computer, and (3) alphabetic information. 

In the case of numerical information, one binary digit of a word is used 
for the sign and 40 binary digits are available for numerical data. Numbers 
are stored in the memory as absolute value and sign. In the arithmetic unit, 
negative numbers may be converted to complementary’ form to keep the 
operational algorithms relatively simple. Thus, negative numbers involved 
in addition, subtraction, and compare are complemented upon arrival in the 
arithmetic unit. In the multiplication, extract, input, and output commands, 
negative numbers are not complemented. 

Numbers may be represented in many different ways. For example, a 
word may represent a signed-binary number lying somewhere between — 2*° 
and +2*, or the binary point may be ahead of the most significant digit in 
which case a word lies in the range —1 to +1; the built-in arithmetic opera- 
tions handle numbers in either of these forms. The word may, on the other 
hand, represent a signed-ten-decimal-digit number where each decimal digit 
is represented as a four-digit-binary number. A floating representation may 
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be used where the first digit represents the sign, the next eight digits repre- 
sent the exponent b, and the next 32 digits represent the significant digits 
of the number in binary form. A floating decimal representation may also 
be used giving numbers of eight significant digits ranging in absolute value 
from about 10-*° to 10*®°. More than one word can be used to represent a 
number to effect much greater precision or range of values. 

Floating operations have been coded (as will be explained later in this 
paper) to provide for the addition, subtraction, multiplication, or division 
of two numbers of the form a-2° (with a and 6 stored in the same address). 
All four floating operations involve about 87 instructions. These floating 
operations are performed in about 3, 3, 6, and 14 milliseconds, respectively. 
Compare these times with 64, 64, and 384 micro-seconds, the times required 
for doing ordinary binary addition, subtraction, and multiplication. 

Instructions are subclassified into three classes. These three classes are 
command words, control words, and code words. 

The command words, commonly called commands, are explicitly ‘“‘under- 
stood”’ and “‘obeyed”’ by the computer. A command causes the computer to 
perform a specific operation and gives the necessary information about where 
to get the needed data and what to do with the results. At present, the 
command system used by the Institute’s computer consists of a set of thir- 
teen commands.’ Eight of the thirteen are what might be termed basic com- 
mands; the other five are variations of these eight. Such a command system 
is known as a four-five address code, with four addresses generally in the 
command word and the fifth address in a control counter. The function, or 
operation, of a particular command is denoted by F, the four addresses of 
the command by a, 8, y, and 4, and the fifth address by e which normally 
determines the address of the next command to be obeyed. 

- The size of the memory determines the length of each address which, in 
the case of a 512 word memory, is nine binary digits. Thus, 36 binary digits 
are used to denote the four addresses, a, 8, y, and 6. Four digits are used to 
denote F; of these four, three are used to define the eight basic command 
words, and one is used to denote modifications of the eight. There is one 
spare digit in the command word. 

The thirteen operations, or functions, of the Institute’s computer, as 
well as the meanings of the addresses for the various command words, are 
given in Table 1. 

In order to change automatically the course of operation in the calculator 
when certain bounds have been reached,‘ there are conditional and uncondi- 
tional transfers of control commands. A conditional transfer is accomplished 
with compare commands. An unconditional transfer is accomplished with 
certain special commands (A;, S;, and M;), wherein the fourth address, 6, 
determines the next command. 

In the table, special compare might well be called absolute compare in 
that the absolute values of the numbers are compared. Since the result of 
the subtraction in the compare operation is put back into the memory, this 
command can be used to obtain the absolute value of a number in one opera- 
tion by comparing the number with zero. The compare command can also 
be used for an operation called tally, as follows: Assume that it is desired to 
repeat a routine fifty times stored in memory locations 31 to 37, inclusive. 
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TABLE 1 
Meanings of the Addresses for the Commands 
Command a B Y 5 F 
Add Address of Address of Address of Address of Next A 
Augend Addend Sum Command if 
Overflow 
Special Address of Address of Address of Address of Next Ai 
Add Augend Addend Sum Command 
Subtract Address of Address of Address of Address of Next S 
Minuend Subtrahend Difference Command if 
Overflow 
Special Address of Address of Address of Address of Next Si 
Subtract Minuend Subtrahend Difference Command 
Multiply Address of Address of Address of M 
Multiplier Multiplicand Product 
Rounded Off 
Special Address of Address of Address of Address of Next Mi 
Multiply Multiplier Multiplicand Product Command 
Rounded Off 
Product Address of Address of Address of Address of Least 
Multiplier Multiplicand Most Sig- Significant Part 
nificant Part of Product 
of Product 
Compare Address of Address of Address of Address of Next c 
Minuend Subtrahend Difference Command if Dif- 
ference is Non- 
negative 
Special Address of Address of Address of Address of Next Ci 
Compare Minuend Subtrahend Difference Command if Dif- 
of Absolute ference of Abso- 
Values lute Values is 
Non-negative 
Extract Address of Address of Address of If Second Digit of E 
Extractor Extractee Extracted 8i ({—shitt Left 
(Determines and Shifted 'S\ 1—Shift Right 
Digits to be Result Other Digits Tell 
Extracted) Number of Places 
to Shift 
Input Address of Drum Address Selects Input I 
Incoming if used Device 
Information 
Special (Incoming Selects Input I; 
Input Information Device 
to 
Address e) 
Output Address of Drum Address Selects Output oO 
. Outgoing if used Device 
Information 


At the beginning of the routine, let the number ‘‘49”’ be placed in address 10. 
Suppose that address 1 stores the number “1.” Place the command 


10, 1, 10, 31, C 


in address 38. The first time this command immediately following the rou- 
tine is obeyed, the number in address 10 will be reduced to “48,’’ and the 
next command will come from address 31. Each time the routine is repeated, 
the command in address 38 will be executed, with a reduction of the number 
in address 10 by one. After the routine has been performed 49 times, the 
number in address 10 will be zero; the 50th time the difference formed by 
the compare order will be. negative, and the next com:nand will come from 
address 39 (the value of ¢). Thus, the desired routine will be repeated exactly 
fifty times. 
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In the case of the normal addition and subtraction commands, overflow 
is automatically detected. An extra digit is provided in the arithmetic unit 
in the most significant end of the A and M® registers so that for normal addi- 
tion and subtraction commands a ‘“‘1”’ in this position will cause the next 
command to come from address 6 instead of address «. In the case of the 
compare command, the proper result of the subtraction is in the A register, 
but here, as in the addition and subtraction commands, the most significant 
digit (overflow digit) is not put back into the memory. In the compare 
command there is no means of detecting overflow. 

The extract command provides for obtaining the logical product of two 
words’ and shifting the result an arbitrary amount; it may also be used to 
delete arbitrary parts of a word. Its primary purpose is to assist automati- 
cally in fabricating new commands during the computation and to sort out 
the exponent from the significant figures in floating-point operations. 

Special input is the command used to insert information into the com- 
puter when it is first put into operation. A word consisting entirely of zeros 
in the F portion is used to designate special input. Thus, when the memory is 
cleared, every word is a special input command, and, with these commands, 
the destination of the incoming data is determined by e. After « has been 
increased by ‘‘1,” it determines the source of the next command. Once 
started, the calculator will count through its complete high-speed memory 
and read in information from the teletype tape.* After the memory has been 
filled, the e« counter steps to zero, and the calculator obeys the command 
stored in that position. 

The input and output commands may specify the magnetic drum as a 
source or destination of words to be transferred. The drum itself will be 
about eight inches in diameter and two feet long and will hold about 8,192 
words.’ Initially only one word at a time can be transferred between it and 
the high-speed memory. The drum will rotate at 3600 revolutions per minute 
which means that the average access time for a word will be 8 milliseconds. 
As soon as possible after the computer is put into operation, counting facili- 
ties will be added to the control to enable the high-speed memory to be 
operated in synchronism with the drum. This will make it possible to transfer 
a whole vector in one revolution of the drum, thus greatly decreasing the 
average access time. 

In general, the drum will serve as an auxiliary storage for numbers, 
instructions, and function tables. The drum is used instead of extra tape 
units because it offers better accessibility to information and requires no 
manual handling. The drum does not hold as much information as a mag- 
netic tape unit, but its size seems adequate for the purposes mentioned 
above. For greater storage several drums may be operated in synchronism.!° 

The second class of instructions, control words, are not directly obeyed by 
the calculator, nor are they a direct part of the calculation; yet in various 
ways, they control the course of the computation or enter into the arithmetic- 
like operations which are performed upon command words. Control words 
may serve as parameters which determine the number of repetitions of cer- 
tain routines; they may be the bounds used to stop certain computational 
processes; or they may serve as factors in the logical products or extraction 
operations. 

The third class of instructions, code words," specify (usually in one word) 
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a whole sequence of procedures for the computer to follow. Thus, one may 
specify, for example, a scalar multiplication with only one code word. This 
code word is made up of parameters which specify the common factor (that 
is, which specify the address in the high-speed memory), the location of the 
elements of the vector (say, by specifying the address of the first element and 
the number of terms in the vector), and the location for the result. 

When prepared for the computer, a problem consists of a sequence of 
words called a main routine. This main routine is made up of instructions 
(commands, code words, and control words) together with the numerical 
constants appropriate to the problem. 

The code words in the main routine contain the parameters or addresses 
necessary to call into action other sequences of instructions, usually called 
subroutines. Subroutines (and code words) exist for such procedures as float- 
ing operations, standard iteration formulae, vector operations, integration 
formulae, and so forth, which are frequently used in the course of doing a 
computation. A subroutine may itself contain code words which call into use 
other subroutines. (All of the more frequently-used subroutines will be stored 
on the magnetic drum, thus being comparatively easily accessible.) 

A routine known as the interpretation routine keeps track of the place in 
the main routine, causes segments of the main routine to be brought into 
the high speed memory, inspects each successive instruction in the main 
routine to see if it is a command (in which case it is obeyed) or a code word 
(in which case the interpretation routine extracts an entry” from the code 
word and sends the computer to the appropriate subroutine). A subroutine 
may itself change the code word being considered by the interpretation 
routine and, thus, cause the computer to carry out other subordinate ac- 
tivities before going on with the main routine. 

Once the appropriate subroutines are established, the task of coding a 
complex problem is very much simpler. For example, the problem of solving 
systems of simultaneous linear equations of orders up to 125 by the elimina- 
tion method has been worked out with the five main subroutines listed in 
Table 2. The arithmetic operations performed by these subroutines are of 


TABLE 2 


Routines for Solving Simultaneous Linear Equations 


Time for 
No. of | Execution 
Name Code Instructions in ysecs. Purpose 
Vector Input a, y, n, VI 8 (+ 76)8000 m words transferred from ad- 


dresses y + i on drum to ad- 

dresses a +4 of high-speed 

memory, i = 0,1, ---,# — 1. 

Vector Output a, y, n, VO 9 (n+ 4)8000 Converse of VI, that i is, transfer 
a! high-speed memory to 


Floating a, B,y,A 88 3,000 Adds, subtracts, multiplies, or 
Operations a, B, 7,8 3,000 divides the "operands and 
a, B, y, M 6,000 puts result in +. 
a, B, y, D 14,000 
Vector Constant a,8,y,",VC 10 6500n +448 Multiplies word in 6 by words 
Product in a +4; answers go into 
y¥+i,i=0, 1, -->,8— 1. 
Vector a, B,y,m,VS 10 3500n +448  Subtracts words in B + i from 
Subtraction those in a + i, answers go into 


y¥+i4,4=0,1,---,n2-—1. 
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the floating binary type, and each such arithmetic operation is in itself a 
subroutine. By using code words the coding for this simultaneous linear 
equation problem is reduced to laying out a sequence of about thirty instruc- 
tions. It is convenient to have all routines and constants, and two rows of 
the matrix stored in the high-speed memory; the figure of 125 is based on a 
512 word memory. The approximate time required to solve a set of sixty 
equations under various conditions is given in Table 3. In fixed-point opera- 


TABLE 3 
Time Required to Solve Sixty Simultaneous Linear Equations 
Fixed Binary Point Floating Binary Point 
Synchronized Non-Synchro- Synchronized Non-Synchro- 
Drum nized Drum Drum nized Drum 
Computing Time 3.5 min. 3.5 min. 19 min. 19 min. 
Transfer Time 1.0 min. 30.0 min. 1 min. 30 min. 
(to and from 
the drum) 
Totals 4.5 min. 33.5 min. 20 min. 49 min. 


tion, a division routine replaces the floating routines of Table 2. The largest 
pivot may be used in each reduction, and scale factors may have to be 


introduced. Table 2, Table 3, and some of the routines were worked out by | 


RosELYN Lipxis of the Machine Development Unit at the Institute for 
Numerical Analysis. 
H. D. HusKEy 


Institute for Numerical Analysis 
Univ. of California, Los Angeles 


1 The Institute is one of four sections of the National Applied Mathematics Laboratories 
of the National Bureau of Standards. It is located on the campus of the University of 
California at Los Angeles. The computing machine discussed in this paper is financed by the 
Air Materiel Command of the United States Air Force. 

*—N is converted to 2 — N or 2° — N, depending upon whether the size of the 
memory is 512 or 1024 words, respectively. 

ese commands are a variation of a set pro by E. F. Moore while he was 
—— in the National Applied Mathematics Laboratories of the National Bureau of 
tan 

‘For example, in the integration of the exterior ballistic equation, the procedure must 
change when a shell has completed its flight; or in a square root iteration, the procedure 
changes when two successive iterants are sufficiently close together. 

* Addition, subtraction, and compare may produce results which exceed the capacity 
of the memory cells. 

6 This extra digit in M is needed for the complement process. 

7 Or it may be used for deleting arbitrary parts of a word. 

8 At a later date, magnetic tape may also be used for inserting information into the 
computer. 

i ® Professor P. Morton of the University of California at Berkeley is constructing such 
a drum. 

10 Professor F. C. WILLIaMs of Manchester University, England, has operated a drum in 
synchronism with a master oscillator. 

11 Code words have been variously termed abbreviated code instructions, quasi-com- 
mands, shorthand commands, abbreviated commands, and coded commands. The term 
“code word” has been selected in preference to these other terms by the author to dis- 
tinguish more clearly this type of instruction from the explicit commands. 

# An entry is the address of the command in the subroutine which should be obeyed first. 

* This table is based on storing three rows of the matrix in the high-speed memory. For 
ap a Cons a 4 sens eighty equations, only two rows can be stored, and input-output 

y fifty per cent. 
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DISCUSSIONS 


Statistical Treatment of Values of First 2,000 Decimal Digits 
of e and of x Calculated on the ENIAC 


The first 2,000 decimal digits of e and of x were calculated on the ENIAC 
by Mr. G. REITWIESNER and several members of the ENIAC Branch of the 
Ballistic Research Laboratories at Aberdeen, Maryland (MTAC, v. 4, p. 
11-15). A statistical survey of this material has failed to disclose any sig- 
nificant deviations from randomness for x, but it has indicated quite serious 
ones for e. 

Let D,‘ be the number of digits 7 (where 7 = 0, 1, ---, 9) among the first 
n digits of e or of x. The count begins with the first digit left of the decimal 
point. If these digits were equidistributed, independent random variables, 
then the expectation value of each D,‘ (with m fixed and i = 0, 1, ---, 9) 
would be 2/10, and the x? would be 


. Ld n\? /n 
= 3 = Se mm —. 
aaa (v ’ is) / 10 


The system of the D,*’s (wherei = 0,1, ---, 9) has 9 degrees of freedom. 
Therefore let 


p = P®(a) = ra [eas 


be the cumulative distribution function of a = x? for k degrees of freedom. 
Then 
bn = P®(a,) 


is a quantity which would be equidistributed in the interval [0, 1], if the 
underlying digits were equidistributed independent random variables. 
Consider m = 2000. In this case, the D,*’s for e are 


(1) 196, 190, 208, 202, 201, 197, 204, 198, 202, 202. 
Hence a, = x,” = 1.11 and p, = .0008. The D,,"s for x are 
(2) 182, 212, 207, 189, 195, 205, 200, 197, 202, 211. 


Hence a, = x,” = 4.11 and p, = .096. 
The e-value of p2ooo is thus very conspicuous; it has a significance level of 
about 1:1250. The x-value of p2000 is hardly conspicuous; it has a significance 
level of about 1:10. 

The relevant fact about the distribution (1) appears upon direct inspec- 
tion. The values lie too close to their expectation value, 200. Indeed their 
absolute deviations from it are 


4, 10, 8, 2, 1, 3, 4, 2, 2, 2, 


and hence their mean-square deviation is 22.2 = 4.71", whereas in the ran- 
dom case the expectation value is 180 = 13.4*. 

In order to see how this peculiar phenomenon develops as m increases to 
2000, a, = xn? and p, of e have been determined from D,‘ for the following 
smaller values of 
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n a, = Xn? Pn 

500 6.72 33 
1000 4.82 BR | 
1100 5.93 25 
1200 i 4.03 093 
1300 3.83 .080 
1400 4.74 .145 
1500 3.69 .070 
1600 2.47 .019 
1700 3.22 .046 
1800 2.85 031 
1900 2.22 013 
2000 1.11 .0008 


These numbers show that the abnormally low value of p, which is so 
conspicuous at » = 2000 does not develop gradually, but makes its appear- 
ance quite suddenly around m = 1900. Up to that point, ~, oscillates con- 
siderably and has a decreasing trend, but at m = 2000 there is a sudden dip 
of quite extraordinary proportions. 

Thus something number-theoretically significant may be occurring at 
about ” = 2000. A calculation of more digits of e would therefore seem to 
be indicated. A conversion to a simpler base than 10, say 2, may also disclose 
some interesting facts. 

We wish to thank Miss Homé McALLIsTER of the ENIAC Branch of the 
Ballistic Research Laboratories for sorting the digital material on which 
the above analyses are based, and Professor J. W. Tukey, of Princeton 
University, for discussions of the subject. 

Since the above was written (November 9, 1949), the ENIAC Branch of 
the Ballistic Research Laboratory very obligingly followed our suggestion 
and calculated the following 500 additional digits! of e. These should replace 
the last 10 digits of the value of e given in MTAC, v. 4, p. 15. 


55990 06737 64829 22443 75287 18462 45780 36192 98197 13991 
47564 48826 26039 03381 44182 32625 15097 48279 87779 96437 
30899 70388 86778 22713 83605 77297 88241 25611 90717 66394 
65070 63304 52795 46618 55096 66618 56647 09711 34447 40160 
70462 62156 80717 48187 78443 71436 98821 85596 70959 10259 
68620 02353 71858 87485 69652 20005 03117 34392 07321 13908 
03293 63447 97273 55955 27734 90717 83793 42163 70120 50054 
51326 38354 40001 86323 99149 07054 79778 05669 78533 58048 
96690 62951 19432 47309 95876 55236 81285 90413 83241 16072 
26029 98330 53537 08761 38939 63917 79574 54016 13722 36188 


This makes it possible to extend the table of a, = x,’ and ~, up to 
n = 2500 


n a, = Xn? Pn 
2100 1.94 .0075 
2200 2.02 .0088 
2300 1.65 .0041 
2400 1.70 .0046 


1.90 
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Thus the values of p, for 2100 < m < 2500 are still significantly low but 
higher than the value of p, at » = 2000. 

Note that the general size and trend of p,, as well as its sudden deviation 
at m = 2000, indicate a non random character in the digits of e. 

More detailed investigations are in progress and will be reported later. 


Los Alamos Scientific Laboratory N. C. METROPOLIS 
Ballistic Research Laboratories G. REITWIESNER 
Institute for Advanced Study J. von NEUMANN 


Princeton, N. J. 


1 Both e and 1/e were computed somewhat beyond 2500 D and the results checked by 
actual multiplication. 


Notes on Numerical Analysis—2 
Note on. the Condition of Matrices 


1. The object of this note is to establish the following theorem. 


THEOREM. Let A be a real n X n non-singular matrix and A’ be its transpose. 
Then AA’ is more “‘ill-conditioned”’ than A. 


This theorem confirms an opinion expressed by Dr. L. Fox! based on his 
practical experience. The term ‘condition of a matrix’’ has been used rather 
vaguely for a long time. The most common measure of the condition of a 
matrix has been the size of its determinant, ill-conditioned matrices being 
those with a “small” determinant. With this interpretation imposed, the 
theorem is clearly correct. More adequate measures of the condition of a 
matrix have been proposed recently by JoHN voN NEUMANN & H. H. 
GOLDsTINE? and by A. M. TurinG.* Their definitions concern all matrices, 
not just the ill-conditioned ones, characterized by very large condition 
numbers. The following two of these definitions will form a basis for the 
proof of the above-mentioned theorem: 


The P-condition number is |Amax|/|Amin|, Where Amex and Amin are the 
characteristic roots of largest and smallest modulus.* 
The N-condition number is N(A)N(A-")/n, where* 


N(A) = (2 aix*)*. 


2. Proof of the theorem in the P case: 

Let A; be the characteristic roots of A and yu; those of AA’ (which are in 
general distinct from the squares of the absolute values of \;). E. T. BRowNE* 
has shown that 7 

min < Ais < pmax. 
From this it follows that 


Amex} < | Amex 


Mmax 
1<¢ — 
= ” Mmin . 














Am in Am in 


which implies the required result. 

3. Proof of the theorem in the N case: 

It is known that N(A) is the square root of the trace of AA’ and therefore 
equal to (> u,)?. The numbers y; are all positive since AA’ is symmetric 
atid positive definite. Since the characteristic roots of A’A and AA’ are the 
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same and since the characteristic roots of the inverse of a matrix are the 
reciprocals of those of the original matrix, it follows that 


N(A~) = (trA-(A-)')! = (tr(A’A))! = (Ew). 
The N-condition number of A is therefore 


1 
rm (XD wi) (Xo wi)! 
In a similar way it can be shown that the N-condition number of AA’ is 


1 
5 (we) wi). 
The theorem follows from the inequality 


Lue Dw? > Dw Da, 
which is in fact true for all real and positive numbers. (It is, indeed, true 
when the first power on the right is replaced by an arbitrary power r and 


the second power on the left by a power s > r.) The proof of the inequality 
is as follows: 


Det Lat — Daw D a 


=n+ Dw? ur *—n— > ens 
ify ix7 


D (utus? + wea) — 2X (wins 1+ pgs) 


i<i 


=D { (ems + we uae 1 + py! — 1) — 2} 2 0 


i<j 


since 
ett oet 22, and py tt wut —1 2 1. 
There is equality if and only if 


ma He°** > 
OLGA TAUSSKY 
NBSMDL 
e an a course of lectures given at the British Admiralty by himself and D. H. SADLER 
in ’ 
* J. voN NEUMANN & H. H. Gotpstine, “Numerical inverting of matrices of high order,” 
Amer. Math. Soc., Bull., v. 53, 1947, p. 1021-1099. (These authors consider symmetric 
matrices only, but it is reasonable to apply the definition to the general case. 


) 
7A. M. TurIne, ‘ a errors in matrix processes,” Quart. Jn. Mech. Appl. 
Math., v. 1, 1948, p. 287-308. 


1E. T. BROWNE, “The characteristic equation of a matrix,’’ Amer. Math. Soc., Bull., 
v. 34, 1928, p. 363-368. 


BIBLIOGRAPHY Z—XI 


1. E. G. ANDREws, “The Bell Computer, Model VI,’’ Electrical Engineering, 

v. 68, 1949, p. 751-756, 7 figs., 5 tables. 22.2 K 29.5 cm. 

Controlled from remote stations, this new digital computer of the relay 
type reduces punched-tape instructions to a minimum. With novel control 
features similar to those used in recent automatic dial-telephone develop- 
ments, this ‘‘upper-class” computer possesses six “‘intelligence levels.’’ Sub- 
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ordinate levels are capable of solving problems such as complex-number 
multiplication without special guidance. 
Author's summary 


2. W. R. Asusy, “Design for a brain,’”’ Electronic Engineering, v. 20, 1948, 

p. 379-383, figs. 

An ideal “thinking machine’”’ must possess negative feedback—or the 
ability to look after itself cleverly by correcting all deviations from a central, 
optimal state. The author describes the homeostat, an electro-magnetic 
machine which is capable of selecting its own arrangement of feedbacks. The 
designer has merely provided it with plenty of variety so that, if the basic 
conditions are altered, it can adapt itself. The author points out that the 
homeostat is still too larval to be a serviceable synthetic “brain.” 

The merits of the ENIAC, the ACE, and the homegstat are compared 
in playing the game of chess. The latter needs no detailed instructions—but 
rather a method by which it is informed of the occurrence of illegal moves 
and mates. Such a machine, if perfect, could in the opinion of the author 
eventually play with a subtlety and depth of strategy beyond that of its 
designer. 

There are two apparent disadvantages to this machine. First, the machine 
will develop a temperament which will probably be manifested in a form too 
complex for the designer’s understanding. More serious in implication is the 
selfishness of the machine; it will judge the appropriateness of an action by 
how it is affected by the feedback. 


NBSMDL Epitx T. Norris 


3. Isaac L. AUERBACH, J. PRESPER ECKERT, JR., ROBERT F. SHaw, & 
C. BRADFORD SHEPPARD, “Mercury delay-line memory using a pulse rate 
of several megacycles,” I.R.E., Proc., v. 37, 1949, p. 855-861, figs. 


A mercury delay line memory system for electronic computers, capable 
of operating at pulse repetition rates of several megacycles per second, has 
been developed. The high repetition rate results in a saving in space and a 
reduction in access time. 

Numerous improvements in techniques have made the high repetition 
rate possible. The use of the pulse envelope system of representing data has 
effectively doubled the possible pulse rate; the use of crystal gating circuits 
has made possible the control of signals at high pulse rates; and a multi- 
channel memory using a single pool of mercury has simplified the mechanical 
construction, reduced the size, and made temperature control much easier. 

The memory system described makes possible a significant increase in 
the over-all speed of an electronic computer. 

Author's summary 


4. Kay Howarp Barney, “The binary quantizer,” Electrical Engineering, 

v. 68, 1949, p. 962-967, figs. 

The binary quantizer is a new device for translating a time-varying 
voltage in an analogue computer into a binary number—the size of each 
being directly proportional to the instantaneous value of the quantity being 
measured. This makes possible the use of digital techniques in an analogue 
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computer. Although primarily designed for computing work, this apparatus 
may find use where quantization of a continuously varying function is 
desired, as in pulse code modulators, automatic metering devices, and re- 
cording machines employing typed or printed numbers. 


5. WARREN H. B.Iss, “Electronic digital counters,”’ Electrical Engineering, 
v. 68, 1949, p. 309-314, figs. 


Recent advancements in all branches of science have created demands 
for high-speed counting devices. The development of the electronic digital 
counters, which are used for many purposes in the large computers, has made 
such fast computations possible. These counters, which use the binary sys- 
tem, are described in detail in this article. 

Author's summary 


6. D. R. HARTREE, Calculating Instruments and Machines, Univ. of Illinois 
Press, Urbana, IIl., 1949, 138 + ix pages, 68 figs., $4.50. 17 K 25.3 cm. 
Professor Hartree has written for the slightly mathematical reader an 

able summary of existing and projected computers and of some of the mathe- 

matical problems which can be solved with their help. His interpretation of 
the subject is authoritative, as he has personally contributed to many of 
the phases of the field of computing machinery. Everyone has been intrigued 
by the Hartree Differential Analyser built of Meccano parts at the University 
of Manchester in England. More recently Professor Hartree has solved 
problems on several of the large-scale digital machines, so that he speaks 
from personal knowledge and experience when he discusses these computers. 

In some respects the book has a distinctly English flavor; the terms 
“instrument” and “‘machine’”’ are used to denote ‘‘analogue”’ and ‘“‘digital” 
computers, respectively. In addition, there are numerous references to de- 
velopments in England which may be less familiar to American readers than 
the parallel developments in this country. 

After an introductory chapter discussing the distinction between instru- 
ments and machines, the author devotes two chapters to the various differ- 
ential analyzers and to their uses in solving ordinary and partial differential 
equations. Chapter 4 concludes the section on instruments with a rapid 
sketch of devices for solving systems of linear algebraic equations, for finding 
the roots of polynomials, and for integration. (It was noted that no mention 
was made of the TrAvis-Hart electrical root-finder.) 

The section on digital machines opens with a much needed chapter on 
terminology. A chapter is devoted to BABBAGE’s astonishing Analytical 
Engine of one hundred years ago, and another to existing computers— 
mechanico-electric, relay, and electronic. Chapter 8 contains schematics 
and discussions of machines now under construction or in the planning stage. 

The book closes with an especially valuable chapter on the methods of 
numerical analysis, emphasizing the ‘‘machine’s-eye view’’ of these methods. 
The author mentions iterative methods and the solution of algebraic, ordi- 
nary and partial differential equations, with useful hints based on his own 
experience. 

Our thanks are due to Professor Hartree for a pleasing and a useful book. 

G. R. Stipitz 


393 South Prospect St. 
Burlington, Vermont 
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7. Harry D. Huskey, “The status of high-speed digital computing sys- 
tems,” Mechanical Engineering, v. 70, 1948, p. 975-977, bibl. 

A brief outline of the history of various computing devices precedes a 
discussion of technical aspects of many of the current large-scale computers 
being built all over the world. The author discusses the réle of computing in 
science emphasizing the continuing necessity for hand computers and desk 
machines as well as the high-speed computers. He goes on to point out that 
the current computers can substitute as a universal model for a large class 
of model experiments. 


Epitu T. Norris 
NBSMDL 


8. T. KiLBurn, “The University of Manchester universal high-speed digital 
computing machine,” Nature, v. 164, 22 Oct. 1949, p. 684-687. 
The author describes the binary high-speed electronic computer being 
developed under the direction of Professor F. C. WILLIaMs, with the active 
assistance of the Telecommunications Research Establishment, Malvern. 


9. WaRREN S. McCuLtocu, “The brain as a computing machine,” Elect ical 

Engineering, v. 68, 1949, p. 492-497, bibl. 

The author, a medical doctor, employs electrical engineering terminology 
to show how the brain may be likened to a digital computing machine con- 
sisting of ten billion relays called neurons. To carry the analogy further, the 
performance of the brain is governed by inverse feedback; subsidiary net- 
works secure invariants, or ideas; predictive filters enable us to move toward 
the place where the object will be when we get there; and complicated servo- 
mechanisms enable us to act with facility and precision. Disorders of function 
are explained in terms of damage to the structure, improper voltage of the 
relays, and parasitic oscillations. 

Author's summary 


10. J. Howarp Parsons, ‘Electronic classifying, cataloging, and counting 
systems,” I.R.E., Proc., v. 37, 1949, p. 564-568, figs. 

The determination of the distribution of a series of physical events ac- 
cording to magnitude is important in the study of the associated physical 
laws. Previous methods of determining this magnitude distribution were 
slow and cumbersome. Three new electronic systems which operate on events 
at a very high rate have been developed by the author while at Oak Ridge 
National Laboratory, and these are described. 

The analyzers can be used to determine the magnitude distribution of 
any series of physical events, if the characteristic under observation can be 
translated into a proportional voltage pulse. Two applications are discussed, 
and the advantages of the analyzers over other systems are shown. 

Author's summary 


11. T. PEaRceEy, “Modern trends in machine computation,” Australian 
Journal of Science, Supplement, v. 10, Feb. 21, 1948, 20 p., 18 figs. 
The object of this paper is “‘to indicate the general structure and organ- 
ization of high-speed computing machines and the avenues which have been 
opened up for elaborate high-speed calculations and to suggest the advan- 
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tages and disadvantages of the various systems available.’’ The author 
briefly discusses the requirements for an ideal high-speed computer, the basic 
organic structure, number representation systems, machine components, 
and simultaneous versus parallel operation. He briefly outlines computer 
developments which were then current and predicts that in the future these 
machines may possibly tackle the symbolic procedures of pure mathematics 
provided the premises are supplied to the machine. The uses of these ma- 
chines are not limited to the solution of mathematical problems; for ex- 
ample, they could be used to set up efficient, speedy filing systems in 
industry. 
Epita T. Norris 

NBSMDL 


12. Tatsujtro Suimizu & Yotcut Katayama, “Solution of non-linear equa- 
tions by punched card methods,” Math. Japonicae, v. 1, 1948, p. 92-97. 
This paper discusses first the solution of simultaneous algebraic equations 

with particular reference to finding the initial approximation to the roots of 

the equation using punched card techniques. After the initial difficulty of 
obtaining the first approximation has been surmounted, closer approxima- 
tions can easily be found using the NEWTON or other iterative techniques. 

The author specifically applies his method to the solution of an ordinary 

differential equation by reducing it to the solution of a system of simul- 

taneous equations. 


Epita T.:- Norris 
NBSMDL 


13. M. V. WiLKEs, ‘Electronic calculating-machine development in Cam- 
bridge,”” Nature, v. 164, Oct. 1, 1949, p. 557-561, illustrs. 


A very brief description of the EDSAC is presented, and the solution of 
a typical problem on the machine, i.e., the evaluation of Arry’s integral is 
discussed. 


14. M. V. WiLkEs, “Programme design for a high-speed automatic calcu- 
lating machine,” Jn. Sci. Inst. and Phys. in Industry, v. 26, 1949, p. 
217-220. 

Problems intended to be solved with the aid of a high-speed digital 
calculating machine must first be reduced to a series of arithmetical opera- 
tions. These, together with various auxiliary operations required for such 
purposes as keeping count of the number of times a particular routine has 
been repeated, must then be expressed in the code appropriate to the ma- 
chine. The paper contains an account of some of the details of this process, 
with special reference to the EDSAC, a high-speed electronic machine in the 
University Mathematical Laboratory at Cambridge. Simple examples are 
given. These have been specially designed to illustrate the use of conditional 
orders and the way in which arithmetical operations may be performed 
on orders. 


Author's summary 
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15. M. V. WiLkes, “Progress in high-speed calculating machine design,” 

Nature, v. 164, 27 Aug. 1949, p. 341-345. 

An outline is presented of computer developments (chiefly those in 
England) which were discussed at a conference on high-speed automatic 
digital calculating machines held in the University Mathematical Labora- 
tory, Cambridge, on June 22-25, 1949, to mark the completion of the 
EDSAC. [For an account of this conference see MTAC, v. 4, p. 51-53.] 


16. R. WiLson WILLIAMs, “‘A survey of some recent advances in computing 
devices,” Science Progress, v. 37, 1949, p. 42-52, bibl. 


The paper presents a concise survey of recent trends in the design and 
use of computing devices. Emphasis is laid on those developments involving 
electronic and electro-mechanical principles. 


NEws 


Institute for Numerical Analysis.—A course on electronic computing machines, with 
special emphasis on the NBS INA Computer, was initiated October 28, 1949, at the 
Institute for Numerical Analysis. This course, which meets once a week under the direction 
of Dr. H. D. Husky, is for the purpose of acquainting interested persons with the possi- 
bilities and limitations of high-speed digital computers of this type. The course is open 
to personnel of other governmental and industrial agencies as well as to UCLA and INA 
staff members. 


Institute of Mathematical Statistics.—The forty-first meeting of the Institute of 
Mathematical Statistics and the twelfth annual meeting was held in New York City on 
December 27-30, 1949, in conjunction with the meeting of the American Statistical Associa- 
tion, the American Association for the Advancement of Science, the American Mathematical 
Society, the Econometric Society, the Psychometric Society, the Mathematical Association 
of America, the Association for Computing Machinery, and the American Psychological 
Association. 

On Friday afternoon, December 30, 1949, a special session of the meeting was devoted 
to computation at which the following papers were presented: 

“Idiosyncrasies of automatically-sequenced digital computing machines’ by Ipa 
Ropes, NBS. 

“Problem solving on large-scale automatic calculating machines” by W. D. Woo, 
Harvard University. 

“A statistical application of the UNIVAC” by Jonn Maucaty, Eckert-Mauchly 
Computer Corporation. 

Discussion by JAMES McPHERSON, Bureau of the Census, and Emm ScHELL, Office of 
the Air Comptroller. 

In the first paper, the speaker stressed the fact that the so-called “‘brain’’ machines are 
in reality morons whose every move must be directed. Hence, problems to be put on these 
machines must be programmed in minute detail; the programmer must foresee every diffi- 
culty which might arise in the problem solution. Dr. Woo described the design features of 
the Mark III calculator, supplementing his discussion with slides. In the UNIVAC discus- 
sion, Dr. Mauchly brought out the fact that the UNIVAC was designed primarily for sta- 
tistical applications; in this machine the emphasis is not on speed for its own sake but rather 
on rapid input and output speeds needed for Bureau of the Census problems. A detailed 
discussion of random number generation on the BINAC was also given. During the discus- 
sion period following these talks, typical problems to be solved on the UNIVAC were listed, 
and the time-saving features of this machine were emphasized. 
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International Business Machines Corporation.—A Seminar on Scientific Computation 
was held at Endicott, New York, Nov. 16-18, 1949 by the International Business Machines 
Corporation. The program for the seminar was as follows: 

“Wednesday Morning, November 16, LEon Brittourn, IBM, Chairman 

“The dynamics of nuclear fission” by Davip L. Hitt, Vanderbilt University. 

“Monte Carlo calculations” by WrLL1AM Woopsury, Northrop Aviation Co. 

“Modifications of the Monte Carlo Method” by Herman Kaun, The RAND Cor- 
poration. 

Wednesday Afternoon, November 16, Frank Hoyt, Argonne National Laboratory, 

Chairman . 

“Analyzing exponential decay curves” by ALSTON HOUSEHOLDER, Oak Ridge National 
Laboratory. 

“Eigenvalue problems related to the diffusion equation” by DonaLp A. FLANDERS, 
Argonne National Laboratory and GEorGE SHORTLEY, Operations Research Office, Depart- 
ment of the Army, Johns Hopkins University. 

Demonstration of the IBM Card-Programmed Electronic Calculator. 

Thursday Morning, November 17, VERNER SCHOMAKER, California Institute of Technology, 

Chairman. 

“Stochastic methods in statistical and quantum mechanics” by GiLBERT KiNG, A. D. 
Little Company. 

“Calculations of resonance energies” by GEORGE KIMBALL, Columbia University. 

“Cam design calculations on the card-programmed electronic calculator” by E. A. 
Barser, IBM. 

Thursday Afternoon, November 17, ARTHUR RosE, Pennsylvania State College, Chairman 
“Distillation theory” by JoHn Bowman, Mellon Institute for Industrial Research. 
“Calculation of multiple component systems” by Stuart R. BRINKLEy, U. S. Bureau 

of Mines. 

Friday Morning, November 18, A. H. Taus, University of Illinois, Chairman 
“The parabolic equation” by L. H. Tomas, Watson Laboratory. 

“Solutions of the wave equation” by Paut HErGEt, Director, Cincinnati Observatory. 

“Sampling methods applied to differential equations” by JoHN Curtiss, National 
Applied Mathematics Laboratories, National Bureau of Standards. 

Friday Afternoon, November 18, W. J. Eckert, Watson Laboratory, Chairman 
“On the distribution of KoLMoGorov’s statistic for finite sample size’”” by Z. W. Brrn- 

BAUM, University of Washington. 

“Specialized matrix calculations on the card-programmed electronic calculator’ by 
H. R. J. Groscu, Watson Laboratory. 


Swedish State Board of Computing Machinery. BARK,-the Swedish relay machine 
[MTAC, v. 4, p. 52-53] built in Stockholm by C. C. R. A. Pat and collaborators under 
the Swedish State Board of Computing Machinery is now completed and. under trial 
running. A description of the machine will be given in a future issue of MTAC. The me- 
chanical differential analyzer, built in Gothenburg by S. EkeLér [MTAC, v. 3, p. 328] and 
collaborators has also been completed, Preliminary planning for an electronic computing 
device is going on under the direction of Ekeléf and Palm. 
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BIBLIOGRAPHY Z-XI 


17. A. E. Carter & D. H. Sapier, “The application of the National 
Accounting Machine to the solution of first-order differential equations,” 
Quart. Jn. Mech. Appl. Math., v. 1, 1948, p. 433-441. [MTAC, v. 3, 
p. 548.] 


Details are given of the solution on the National machine of the first- 
order differential equation y’ = f(x, y) by the MILNE-STEFFENSEN method. 
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The formula used is Yn — yas = (3g + 282 + yyo'g)n-2, where g(x, y) 
= $hf(x, 9). 

The machine is used to build up (3g + 2é?g),_2 from initial values and 
5‘g up to 5‘g,_3, to build up an approximate y, from yn_«, (3g + 28g),-2, and 
and extrapolated 4%5‘g,-2, and to difference the approximate g, to 5‘g,_». The 
operator has to provide the extrapolated ,46‘g,_2, to obtain the approximate 
ga, and to correct y, and *g,_» by successive approximation. He must further 
enter these values in the machine at the appropriate times. 

The machine prints on a line, m, y,(approx.), ga—1, ga(approx.), d*g,_s 
(approx.), S*gn2, 5*gn—s/2, Fgn—1, 5En—1/2, Ya-s, and yy5*g,1(extrapolated). 
The second and last of these are corrected by hand. 

The next correction term’ of the series, —z?,5°g, can be applied to y by 
hand if necessary. 

Anyone who wishes to use the Milne-Steffensen method should be 
warned that it introduces logarithmic rates of increase of error in the series 
of y, additional to that, h df/dy, inherent in the differential equation. In this 
x. h of 19 af 
3° + ay *  - Sy: 
Thus, it may be unwise to use the method to integrate in the direction for 
which df/dy is negative. -_ oe 


case, the additional logarithmic rates are + 


Watson Scientific Computing Laboratory 
612 West 116th Street 
New York 27, N. Y. 


18. A. C. Cook & F. J. Macrnniss, “More differential analyzer applica- 

tions,”’ Gen. Elec. Rev., v. 52, no. 8, 1949, p. 14-20. 

This article lists a variety of engineering applications of the G.E. differ- 
ential analyzer, which has been described in articles by PETERSON & KUEHNT' 
and by Peterson & Concorp1A.? The applications include electron ballis- 
tics, long distance power transmission, bearing design, frequency changers, 
guided missiles, and control problems. These applications extend those given 
in a previous paper by MAGINNIss.* F.J.M 


1H. A. Peterson & H. P. Kuvennt, “A new differential analyzer,” A.I.E.E. Trans., 
v. 63, 1944, p. 221-228. [MTAC, v. 1, p. 430-431.] 

2H. A. Parmmsost & C. Concorp1a, “Analyzers for use in engineering and scientific 
problems,” Gen. Elec. Rev., v. 48, no. 9, 1945, p. 29-37. [MTAC, v. 2, p. 55.] 

3F. J. Macrnniss, “Differential analyzer applications,” Gen. Elec. Rev., v. 48, no. 5, 1945, 
p. 54-59. [MTAC, v. 1, p. 452-454. ] 


19. A. B. Macnee, “A High Speed Electronic Differential Analyzer, I.R.E. 

Proc., v. 37, 1948, p. 1315-1324. 

A differential analyzer is described which reproduces its solutions sixty 
times a second, thus permitting the results to be exhibited on a cathode ray 
tube. Resetting is accomplished by clamping circuits. Addition and integra- 
tion are performed in the usual fashion (involving feedback amplifiers). The 
multiplier uses a cathode ray tube. An electron beam is bent by electrostatic 
and magnetic methods to yield an effect proportional to the product of two 
quantities. The resultant deflection is cancelled by a feedback controlled 
electrostatic defiection, which measures the product. When the feedback is 
applied to a factor, division is obtained. An arbitrary function is introduced 
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into the machine by a mask on the face of a cathode ray tube. A feedback 
arrangement causes the beam to follow the edge of the mask. 

An error analysis for differential equations with constant coefficients is 
also given. The results of applying this differential analyzer to a variety of 
problems are given including end point boundary problems and the equa- 
tions of MATHIEU, HILL and VAN DER POL. F. J. M 


20. J. C. JAEGER & J. D. CLarKE, ‘“‘A Product Integraph,” Jn. Sci. Instr., 
v. 26, 1949, p. 155-156. 


By using mechanical components obtained from surplus anti-aircraft 
predictors, the authors have built an instrument intermediate in accuracy 
and cost between Meccano constructions and the well-known large differ- 
ential analyzers. Two ball-and-disk integrators permit the evaluation of the 
indefinite integral, 


fF terete. 


The input curves f(x) and g(x) are followed manually; the machine is motor- 
driven in the x direction: 
H. R. J. Groscu 


Watson Scientific Computing Laboratory 
612 116th Street 
New York 27, N. Y. 


21. L. PEREK, ‘‘Nomograms for computing galactic longitude and latitude, 
galactic rectangular coordinates and components of the space velocity,”’ 
Brno, Masarykova Univ., Prirodovéck4 Fakuta, Spisy, no. 299, 1947. 
The first two nomograms constructed for finding the galactic longitude 

and latitude of a star are at present of little practical importance, since 

detailed tables for the conversion of right ascension and declination to 
galactic longitude and latitude have already been published.! 

The second set of seven nomographic charts is designed for the calcula- 
tion of the rectangular components x, y, z of space velocities. The coordinate 
system is so oriented that the z axis points to the galactic North Pole 
(RA 12540™, decl. +28°), the x axis toward the center of the Galaxy (galact. 
long. 327°). The data from which the velocity components are derived are 
the proper motion components ua, ws, the parallax w (or distance 7), the 
radial velocity V, the right ascension a and declination 6 of the star. The 
nomographic solution of this problem is quite accurate enough, as space 
velocities of stars can generally not be determined with an accuracy greater 
than 1%. The procedure, however, is by no means simple; it requires eleven 
nomogram readings of which seven are double constructions. So long as 
extensive tables are not available these nomograms may, nevertheless, 
serve a useful purpose. 

Three of the charts made for the calculation of velocity components 
serve also to obtain the rectangular space coordinates of a star (in the same 
system) from the right ascension, declination, and distance. 

R. J. TRUMPLER 

University of California 

Berkeley, California 


1 J. Ontsson, “Tables for the conversion of equatorial coordinates into galactic coordi- 
nates,” Lund Observatory, Annals, v. 3, 1932. 
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22. RoBERT M. WALKER, “An analogue computer for the solution of linear 
simultaneous equations,” I.R.E., Proc., v. 37, 1949, p. 1467-1473, figs. 
Linear simultaneous equations occur frequently in science and in en- 

gineering. Their solution by numerical methods is straightforward, but the 

amount of work required increases rapidly with the number of unknowns. 

A device is described for the solution of systems of linear simultaneous equa- 

tions with not more than twelve unknowns. It is an electrical analogue com- 

puter which accepts the problem information in digital form from a set of 
punched cards. This facilitates the preparation, checking, and insertion of 
the input data and greatly reduces some of the usual liabilities of an analogue 
device. No special preparation of the problem is required, other than a simple 
one of scaling the coefficients. Solutions of well-determined problems are 
easily and rapidly attained and may be refined to any desired accuracy by 

a simple iteration procedure. 

Author's summary 


NOTES 


112. A CoMMITTEE ON Factor TABLES.—In September, 1946, the Asso- 
ciation Frangaise pour l’avancement des Sciences established a committee 
consisting of A. GERARDIN (France), who to our regret could not participate 
in the work because of ill health, L. Potetti (Italy) and the author, for 
the purpose of extending the factor table. The committee was joined later 
by Dr. A. GLopEN (Luxemburg). We agreed that only a table practically 
free from error may have any significant value. Hence, the necessity to check 
the existing manuscript tables against one another. These tables are: 
KULIk’s famous manuscript,' Poletti’s table of the 11-th million, GoLUBEV’s 
table of the 11-th and 12-th millions, R. J. PoRTER’s two tables? of the 11-th 
million. The Carnegie Institution of Washington presented us with two 
microfilms of the 11-th and 12-th millions of Kulik’s manuscript. Poletti’s 
table is in our hands. A request of a photostat of Golubev’s table in the 
Steklov Institute at Moscow finally was denied. Porter’s tables are unfortu- 
nately in symbols quite different from Kulik’s. The committee decided to 
extend the existing printed table to the 11-th and probably the 12-th million. 
The necessity of checking Kulik’s table against Poletti’s made it necessary to 
get large photos of the relevant part of the microfilm. The Lord Mayor of 
Luxemburg presented us with large photos of the second half of the 11-th 
million of the microfilm and the author ordered large photos of the other 
half. Since Kulik used letters instead of numbers, the author is to transcribe 
the photos in the interval from 10’ up to 10’ + 5-10°, Dr. Gloden in the 
interval 10’ + 5-105 up to 10’ + 10°. Every page of the new manuscript 
will be checked against Poletti’s table. The cliecking of the manuscript 
against a third table, which seems necessary is still a problem to the 
committee. 

N. G. W. H. BEEGER 
Nicolaas Witsenkade 10 
Amsterdam Z, Holland 


1See MTAC, v. 2, P 139-140, v. 3, p. 222. 
2 MTAC, v. 1, p. 451. 
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113. SupERsoNIc FLow CatcuLations.—The following is an account 
of work done on the ENIAC on supersonic flow past cone cylinders. 


Numerical solutions for equations for the flow of a compressible gas at 
supersonic speeds past a cone cylinder with attached shock wave at various 
combinations of Mach number and cone semiangle have been computed 
on the ENIAC [MTAC, v. 3, p. 206-207]. The general method of obtaining 
these solutions will be described first. Using characteristic variables, the 
equations of irrotational motion for the supersonic flow past a body of revo- 
lution may be written as follows: 


Aya — (K + R)xa 

Hyg — (K — R)xs 

Hua + (K — R)ta + (a*v/y)xa 
(K — R)ug + Lug + (a*v/y) ye 


where x, y are cylindrical coordinates with the x axis along the axis of the 
projectile, u, v are reduced velocity components in the direction of the x, y 
axes, where u = i/c, v = i/c, and c is the limiting velocity as the velocity 
of sound approaches zero, #, i are the original velocity components, a? = a/c? 
is the similarly reduced velocity of sound: 


H=a@—- 7, K = — wW 
L=a-?, R= a(¢ — a’). 


The surface of the body is assumed to be a stream surface, and with the 
RANKINE-Huconiot shock wave conditions give us the following boundary 
conditions, in which the subscripts 1 refer to conditions ahead of the shock 
and the subscripts 2 refer to conditions immediately behind the shock: 


y = F(x) is the contour of the given body, v = uF’(x), 
vPE(y — 1)/(y + 1)qi + 2q:/(y + 1) — ue] 
= (gi — u2)*Lu2 — (y — 1)/(v + 1)qr] 


dy/dx = (qi — tz)/v2, where dy/dx is the slope of the tangent. 

y = 1.4 is the ratio of the specific heats. 

The values of x, y, u, v were computed at the intersections of all charac- 
teristics a = const., 8 = const. in the. following region. The characteristic 
a = 0 was chosen as the downstream characteristic emanating from the 
shoulder in TaYLOR-MACcCcOLL flow about the cone. Along this characteristic 
the values of x, y, u, v were computed at k intervals A@ terminating at the 
shock wave. 

The characteristic 8 = k emanates from the intersection of a = 0 and 
the shock wave, and the region under discussion is bounded by a = 0, 8 = k, 
and the generating curve of the cylinder. In this region the exact flow is 
irrotational since the shock wave is straight. 

The shoulder of the cone-cylinder is mapped onto the characteristic 
8 = 0, and PRANDTL-MEYER flow is computed along this characteristic, 
using k points. The generating line of the cylinder is mapped onto the line 
B = a — kin the a, 8 plane. In the computation of the Taylor-Maccoll flow 
the partial differential equations were reduced to a set of four simultaneous 
ordinary differential equations, with 


AB = k(v/u)/[(v/u)cone — (v/u)shock]. 
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Likewise the equations were reduced to ordinary differential equations 
for the Prandtl-Meyer flow, with Aa = k(v/u)/(v/u)cone. 

The first computations carried out were for M = 2.1297 and a cone 
cylinder of 20° semi-angle, with k = 2, 3, 4, 5, 6, 8, 10, 12, 15, 16, 20, 24, 30, 
32, 40. Taylor-Maccoll flow was computed with k = 960 and Prandtl- 
Meyer flow was computed with k = 16, 32, 64, 128, 256. The partial deriva- 
tives were replaced by difference quotients of first, second and third order 
approximations, and the results were extrapolated to grid size zero (k = ©). 
By using the extrapolation formula z = 2: + (22 — 2;)/3, where z = x, y, u, 9, 
and the subscript 2 indicates that the & for that computation is twice the k 
for the other, it was found that results could be obtained to a much higher 
accuracy in a short time. 

Accordingly, in the next calculation k = 8 and 4, and the values of x, y, 
u, v were calculated at all intersections of the characteristics as well as the 
values of p/p1, p/p1, T/T; at all surface points of the body, where p, p, T stand 
local pressure, density, and temperature, and the subscript 1 indicates free 
stream values. 

The various combinations of cone-angle @ and free stream Mach number 
M for which computations were made are as follows: 


@= 5° M = 1,3, 1.5, 1.7, 2, 3, 4, 5, 7 
@= 9.5° M = 1.3, 1.5, 1.7, 2, 3, 3.8 
6 = 10° M = 1.3, 1.5, 1.7, 1.72, 2, 3, 4, 5, 7 


@ = 12° M =1.7 

@ = 15° M = 1,3, 1.5, 1.7, 2, 2.13, 2.3, 2.7, 3, 4, 5, 7 
6 = 20° M = 1:5, 1.7, 2, 2.13, 3, 4, 5,7 

@ = 25° M = 155, 1.7, 2, 2.13, 3, 4, 5, 7 

6 = 30° M = 1.7, 2, 2.13, 2.3, 3, 4, 5, 7 


6 = 35° M = 2, 2.13, 3, 4, 5, 7 
@ = 40° M = 3,4,5,7 

@ = 45° M = 3,4, 5,7 

@ = 50° M =4,5,7 


The results are in some places good to 4 significant figures and it is 
believed that they are everywhere good to three significant figures. 
R. F. CLIpPINGER 
Ballistic Research Laboratories 
Aberdeen Proving Ground, Md. 


114. NBSCL TaBies.—Eleven of these tables have appeared in second 
editions. We have already referred in MTAC to two of these: Tables of the 
Bessel Functions Jo(z) and J,(z) for Complex Arguments, 1943 (v. 3, p. 25) 
in 1947; and Tables of the Exponential Function e*, 1939 (v. 3, p .173), in 1947. 
There were also a second printing of Table of Circular and Hyperbolic Tan- 
gents and Cotangents for Radian Arguments, 1939 (v. 3, p. 88), in 1947; a 
second edition of Tables of Circular and Hyperbolic Sines and Cosines for 
Radian Arguments, 1939 (v. 1, p. 45), in 1949; and a second edition of Tables 
of Probability Functions, v. 2, 1942 (v. 1, p. 48), in 1948. 

With the approval of the editor of the Journal of Mathematics and 
Physics, the following seven short tables in the MT series were reissued by 
the Government Printing Office to meet a continuing demand: 
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MT 19, “On the function H(m, a, x) = exp (—ix)F(m + 1 — ia, 2m + 2; 
ix),”’ 1942 (v. 1, p. 156) on 15 June 1949; 

MT 20, “Table of integrals /o7 Jo(t)dt and Jo* Yo(t)dt,” 1943 (v. 1, p. 154), 
10 Sept. 1948; 

MT 21, “Table of Jio(x) = S2” Jo(t)dt/t, and related functions,” 1943 (v. 1, 
p. 155), on 25 Feb. 1949; 

MT 22, “Table of coefficients in numerical integration formulae,” 1943 (v. 1, 
p. 157), on 15 Dec. 1949; 

MT 23, “Table of Fourier coefficients,” 1943 (v. 1, p. 192), on 25 Nov. 1949; 

MT 25, “Seven-point Lagrangian integration formulae,” 1943, on 1 June 
1949; 

MT 27, “Table of coefficients for inverse interpolation with central differ- 
ences,” 1943 (v. 1, p. 126, 359), on 1 June 1949. RCA 


115. On LarGE PRIMES AND FACTORIZATIONS.—The following results of 
extensive calculations have been reported. 

I. The number N = 3-2” + 1 is composite. 
It was found that 


2-1 = 7949 31660 65193 00342 32702 86136 93493 91372 (mod JN). 


If N were a prime, we would have 2”—! = 1 (mod N). 

The number N was suggested as a likely candidate for primality by 
Mr. THOROLD GosseEt of Cambridge, England. Thus another attempt to 
discover a larger prime than 2! — 1 ends in disappointment. 


H. S. UHLER 
206 Spring St. 
Meriden, Conn. 


II. The number 6” + 1 is completely factored. 
6% + 1 = 37-313-2341-629 19466 95217. 


The largest factor was proved to be a prime by a method based on the 
converse of FERMAT’S theorem, as described in MTAC, v. 3, p. 496-497; 


v. 4, p. 54-55. N. G. W. H. BEEGER 
Nicolaas Witsenkade 10 
Amsterdam Z, Holland 


III. (a) The number N = (2® + 1)/3-179 is composite 
(b) 3% + 1 = 2-41-6481-28 24290 05041 
(c) 3% + 1 = 2-5-956353- 1743831169 
(d) 2% + 15 is a prime. 
The result (a) follows from the fact that if y = 2®, 


356 — 3 = 3303 41699 82572 42322 50798 (mod JN). 


If N were a prime this remainder would have been zero. The complete 
factorization of 2° + 1 would be of interest since it would probably lead to 
some new multiply perfect numbers. 

The results (b) and (c) fill in two blank entries in CUNNINGHAM & 
WoOOoDALL, p. 11. The former result was found by the method of MTAC, v. 3, 
p. 96-97. The latter result was found by expressing (3% + 1)/10 as the 
difference of two squares. 











It is 


are 
con: 
kno 


Collé 
Allie 


one 
fam 


105 | 
Wal, 


1934 


“for 
wou 
wou 
vali 
tab] 
fixe 
sin : 
real 


whe 
as 71 


the 


and 





of 


to 


he 
7 ; 








QUERY—QUERIES—REPLIES—CORRIGENDA 125 


The result (d) was another example of the converse of Fermat's theorem. 
It is easy to see that the 15 integers between the two primes 


fi, = 22 —1 and po = 2% + 15 


are composite, being multiples of small primes <13. Hence, p; and p» are 
consecutive primes, and constitute the largest pair of consecutive primes 
known. 


A. FERRIER 
Collége de Cusset 


Allier, France™ 


1A. J. C. CunnincuamM & H. J. Woopa., Factorization of y" + 1, London, 1925. 


QUERY 


34. ARCHIMIDES CATTLE PROBLEM.—Has any attempt been made to use 
one of the modern electronic computing devices to get a solution of the 
famous cattle problem! of Archimedes? 


E. P. ADAMs 
105 Plimpton St. 


Walpole, Mass. 


1 See L. E. Dickson, History of the Theory of Numbers, Washington, 1920 and New York, 
1934, v. 2, p. 342-345. 


QUERIES—REPLIES 


44. TABLES OF sin nx/sin x (Q16, v. 2, p. 61).—A table of this function 
“for large integral values of m, say up to 100, and for values of x in radians” 
would be impossibly large if it were interpolable in x. Millions of entries 
would be necessary to carry the table as far as x = 10 with 8D. Since isolated 
values of this function can be obtained with a little trouble from a good 
table of sin x it is not surprising that there is no such table in print. For 
fixed x, however, there are a number of small tables of what is essentially 
sin mx/sin x. Two examples may be cited, though they do not correspond to 
real values of x. For x = arccos(—i/2) we have 


sin nx/sin x = 7'-*F,, 
where F,, is the n-th term of the Fibonacci series, which is tabulated as far 
as nm = 128 [MTAC, v. 2, p. 343]. For x = arccos (—3i/V2) we have 
sin nx/sin x = (iv2)'-"(2" — 1) 
the values of which can be found easily from a table of powers of 2. 


As is well known, U, = sin nx/sin x isa special kind of Lucas! function 
and can be computed recurrently by the formula 


Uae = 2cosxU, — Uy-1. 


- o. 

1. Lucas, Théorie des Nombres, Paris 1890, p. 319. 

CORRIGENDA 
V. 3, p. 362, 1. —1 for RMT 593 read RMT 592. 
V. 3, p. 554, 1. 3, for 4-17, read 3-16. 
V. 3, p. 559, 1. 4, for 454 read 554. 
V. 3, p. 562, 1. 30, for NicHoLas DE Cusa read NICHOLAS DE Cusa. 
V. 3, p. 563, 1. 6, for fifteenth read thirteenth. 
Nos. 28, 29, cover 3, interchange lines J and I. 








